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ABSTRACT 


Acceleration using centrifugal relays is a recently discovered method for 
the acceleration of spaceborne payloads to high velocity at high thrust. Cen- 
trifugal relays are moving rotors which progressively accelerate reaction mass 
to higher velocities. One important engineering problem consists of accu- 
rately tracking the position of the projectiles and rotors and guiding each 
projectile exactly onto the appropriate guide tracks on each rotor. The top- 
ics of this research are the system kinematics and dynamics and the com- 
puterized guidance system which will allow the projectile to approach each 
rotor with exact timing with respect to the rotor rotation period and with 
very small errors in lateral positions. Kinematics studies include analysis of 
rotor and projectile positions versus time and projectile/rotor interactions. 
Guidance studies include a detailed description of the tracking mechanism 
(interrupt of optical beams) and the aiming mechanism (electromagnetic 
focussing) including the design of electromagnetic deflection coils and the 
switching circuitry. 
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Chapter 1 


INTRODUCTION 


Three billion years ago, organic life began to develop rapidly on spaceship 
Earth. It was powered by energy of nuclear origin transmitted from the 
sun in the form of light. Large amounts of that solar energy were then 
stored in the form of coal, oil, and gas foUowing the death of trees and other 
lmng organisms. More recently, man appeared on the spaceship. Until 
about 200 years ago or so, men lived in small separate groups. Then they 
began to integrate their activities into modern societies making the need 
for energy bigger and bigger until they reached a state where those easy- 
to-get hydrocarbons are almost gone. The days of cheap fossil energy are 
most over [1], Nuclear fission energy cannot be more than a transitional 
scheme since the fuel supply is limited and the question of nuclear waste 
has not been solved. People are now turning to solar energy and thinking 
about going to outer space where it can be easily harvested. This idea first 
appeared in 1968 when Peter Glaser proposed a project for orbital space 
power satellites [2], In the mid-seventies, Gerard O’Neill developed further 
the concept of space industrialization [3,4], and by the iate-seventies, the idea 
had gathered momentum and consideration among industry and government 
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m the United States, Europe, and the former USSR [1]. 

As we desperately search for new sources of energy, space industrializa- 
tion, with heavy production of aluminium, iron, silicon, and other energy 
consuming products, appears as an encouraging way out of the current cri- 
sis [1]. Recent years have seen encouraging progress in studies of the trans- 
port, smelting, and use of nonterrestrial resources for large-scale space in- 
dustrialization. A promising approach to transporting such resources is the 
launch from the moon of unprocessed lunar ore. The ore packages would 
be captured and processed in space [5]. Metallic ore can be obtained not 
only from the moon surface, but also from asteroids and it can be processed 
in high Earth orbit as suggested by O’Neill in 1976 [6]. The products can 
then be used for the construction of space structures or sent back to the sur- 
face of the Earth casted into re-entry bodies [1]. It appears though that the 
utilization of lunar materials will almost certainly be a near term necessity 
while asteroidal resources will probably not be tapped for at least another 
ten years according to the most optimistic estimates [7]. The need for fully 
developed industrial schemes both in space and on the surface of the moon, 
and for sophisticated propulsion systems is then explicit. 

1.1 HISTORICAL CONTEXT 

Current launch vehicle systems such as the Space Shuttle, involve launch 
costs on the order of $1000/kg of payload [8]. Before the large-scale de- 
velopment of space resources and energy can succeed, launch costs must be 
reduced. Launch costs for liquid-fuel rockets cannot be lower than the cost of 
the fuel, which is about $3 /kg for liquid hydrogen [9], This limits the cost of 
launching to a minimum of $50/kg of payload, not considering vehicle devel- 
opment and launch operation costs [8], This cost may be too high for cheap 
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space industrialization. Various schemes have been proposed for low-cost 
launch. L. W. Jones et al. proposed a laser propulsion system [10], where a 
ground-based laser would provide the external power. A space elevator with 
electrical power distributed along its length called “The Orbital Tower” was 
suggested by J. Pearson [11], Y. Artsutanov et al. [12] and H. Moravec [13] 
proposed a rotating tether to pick up large payloads from the ground or upper 
atmosphere. A launch loop or “Orbital Ring,” with rapidly moving internal 

components to support it against gravity was discussed by P. M. Birch [14] 
and K. Lofstrom [15]. 

All these different schemes are attractive theoretically, but they require 
either great improvements in materials or lasers, enormous masses in orbit, 
or large, dynamically stabilized masses overhead [8], A compromise system 
putting most of the energy requirements on the ground and requiring mini- 
mum mass in orbit has been suggested by J. Pearson [8] in the late-eighties. 

nsists of a relatively small electromagnetic gun on the ground with a rel- 
atively low-mass rotating tether in orbit to catch the projectile. This system 
requires high launch accelerations and is then limited to unmanned payloads. 

Another concept that got the attention of many scientists is the princi- 
p of solar sails. They have been the subject of several studies on external 
momentum sources for accelerating objects in space. The momentum is ob- 
tained when photons bounce off a reflector [16], Lasers have been suggested 

33 an alternative ^ the sun [17], but the problem is that photons are ineffi- 
cient for this purpose and lead to low momentum/energy ratio. The photon 
is a poor choice of momentum transfer. An additional problem is the quan- 
tum mechanics limitations. It is difficult to collimate a beam adequately. 

Therefore, solar sails are limited to very low thrust and simple laser-driven 
sails are impractical [18]. 
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To overcome the limitations of solar and laser-driven sails, the use of 
macroscopic objects to transfer momentum from a remote source appeared 
more promising. A kinematic study to this approach has been worked out [19] 
and the requirements for maintaining a collimated stream of projectiles is 
achievable [20]. The unsolved problem was the production of a suitable 
stream of high velocity projectiles. Two methods of launching such projec- 
tiles from extraterrestrial platforms were investigated. The first one is the 
linear induction motor [21, 22]. Unfortunately, the required length of a lin- 
ear induction motor increases with the square of the launch velocity, making 
the high-velocity launcher very cumbersome. The second method consists of 
using centrifugal launchers. 

1.2 CENTRIFUGAL LAUNCHERS 

As an alternative to both chemical and elctromagnetic launchers, centrifu- 
gal projectile launchers were proposed to provide reliable, efficient, compact 
systems that can accelerate projectiles to 2 — 3 km/s with energies of 10 5 to 
10 6 J without chemical explosives [23]. The direct conversion of rotational 
mechanical energy into transational mechanical energy with no intermedi- 
ate stores or transfers, makes the potential efficiency in repetitive operation 
higher than any other system. Otherwise, a homopolar generator/rail gun 
system must convert mechanical energy to electrical, store the electrical en- 
er gy inductively, convert the electrical energy to plasma heat and kinetic 
energy, and finally transfer the plasma kinetic energy to projectile kinetic 
energy [23]. The compactness is limited mainly by the yield strength, which 
is the ultimate limitation to compactness of any solid projectile launcher. 

The use of centrifugal launchers in space had been suggested for lifting 
lunar material resources off the moon in the late-seventies [21, 5], but rotors 
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capable of the required 2.2 km/s launch had not been developed yet, so long 
linear induction launchers were chosen as the primary launch option [24]. 

During the fall and winter of 1976, H. Kolm and G. O’Neill made a device 
that would operate at high acceleration (100 g) necessary for a mass driver 
to be used in space. The model served to educate people on the principle 
of mass driver design and also demonstrated that such acceleration could be 
achieved with fairly simple machinery [25]. 

The performance of centrifugal launchers became more known as further 
studies were conducted and the method is now widely used for refuelling a 
plasma in toroidal magnetic confinement devices such as tokamaks and stel- 
larators/heliotrons [26]. Rotors injecting frozen hydrogen pellets into hot 
plasmas were tested at rotor tip speed of 1 km/s. As the fragile hydrogen 
pellets tend to break up at somewhat lower tip velocities, there seemed to 
be little incentive to develop higher tip speeds theoretically allowable by me- 
chanical stress limitations [24]. The proposed lunar application to centrifugal 
launchers is for launching large sacks of raw dirt and the existing type of high- 
velocity centrifugal launcher is best suited to material refined into small balls. 
Therefore, these centrifugal launchers are not feasible for lunar application 
but they may be useful for propelling spacecraft payloads [24], 

Unfortunately, stress considerations limit simple centrifugal launchers to 
a launch velocity not greatly exceeding the exhaust velocity of chemical rock- 
ets [18]. Recently, however, a method has been proposed for substantially 
augmenting the projectile velocities using centrifugal launchers [24]. 

1-3 CENTRIFUGAL RELAYS METHOD 

The basic idea starts with an existing centrifugal launcher [27] as illustrated 
in Figure - 1 [24]. The mass driver is anchored on a large body and used 
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Figure - 1, Centrifugal launcher. 
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to throw small masses towards a curved track or guide tube on the payload. 
Momentum may be transferred by stopping these masses at the payload, or 
more efficiently, by reflecting them back to the launcher [24]. These two 
options are summarized in Figure - 2 [24]. 

To increase the velocity of the projectiles incident on the payload, a cen- 
trifugal relay is inserted between the launcher and the payload as illustrated 
in Figure - 3 [24]. The requirements on this relay are: its rotational energy 
must be maintained so that it accelerates incoming projectiles to a higher 
speed when passing them on to the payload; and it must be continuously 
accelerated so that its velocity is intermediate between the first launcher and 
the payload. This method can be extended by inserting two or more centrifu- 
gal relays between the primary launcher and the payload [18]. This makes 
this recently discovered method using centrifugal relays for the acceleration 
of spaceborne payload to high velocity at high thrust, very attractive. It 
has unique advantages over chemical rockets (low velocity - high thrust), 

and other advanced methods such as plasma thrusters and solar sails (high 
velocity - low thrust) [18]. 

1.4 OBJECTIVE 

The purpose of this work is to explore the use of centrifugal relays for space- 
craft propulsion in sufficient detail to assess the potential of this method for 
orbit transfers requiring large velocity increments (Ay’s up to 10 km/s), and 
for the capture into Earth orbit of suborbital payloads. The investigation of 
this potentially viable novel spaceborne propulsion concept is a joint project 
between the Aeronautical and Astronautical Engineering department, and 
the Nuclear Engineering department at the University of Illinois. This new 
initiative in advanced propulsion systems is funded by NASA. 
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Figure - 2. Ball-payload interation for 
a) stopping and 
b) reflection. 
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Figure - 3. Centrifugal relay method. 




1.5 DEFINITION OF THE PROBLEM 

For the purpose of this thesis, the goal is to assess the potential of centrifugal 
relays for a specific type of mission. The research can be divided into four 
areas: kinematics, mechanical design, guidance, and systems. We propose to 
investigate two aspects of this project and identify all the requirements for 
the kinematics and the guidance system. We will then identify the concepts 
of the solutions to each problem and we will define the components of each 

part. 

We will start in Chapter 3 by giving an overview of the launch system as 
a whole, followed by a description of the kinematics involved to understand 
how momentum is transferred. We will also discuss design considerations for 
the rotors. We will then give in Chapter 3 a detailed description of the launch 
of one ball. Then, we will add to it the catching part in Chapter l We will 
also consider in this chapter the trajectory of the ball on the catching rotor 
and investigate a delay circuit that wiU insure that the ball flies off this rotor 
in the desired direction (which is, in this case, back to the launcher). After 
this one ball, one launcher, and one relay (payload in this case) problem is 
solved, we can consider the case of launching two or more balls. Then we 
can study the problem of additional tracks on the relay. The purpose would 
be to launch balls in the opposite direction and insure conservation of the 
angular velocity of the rotor. This generalization process of launching n balls 
is described in Chapter 5. We also develop in this chapter a program that 
keeps track of the position of the relay as the balls are launched. It comp.les 
all the data necessary for the visualization of the process, such as the time at 
which a ball is released near the center of the primary launcher, the time at 
which it exits the launcher, the time at which it is caught, sent back to the 
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launcher, and at what speed the relay is moving. Finally, this rotor is to be 
considered as a launcher and the sequence repeats itself. The final step would 
be to write a general code that treats arbitrary numbers of relays taking into 
account the externally imposed gravitational potential appropriate to the 
mission type, but this is not in the scope of this thesis. We are limiting 
ourselves to the complete study of a one relay case. For guidance purposes, 
we will be considering guide barrels as a first step in velocity trimming. 
Electromagnetic focussing has been identified as the aiming mechanism. We 
will be giving in Chapter 6 a detailed description of the tracking mechanism 
(interrupt of optical beams) and the switching circuitry for the deflection coils 
and all the electronics involved in the guidance system. The final step would 
be to consider the system as a whole for optimization of mission profiles. One 
could study the Mars mission. One could estimate the performance and the 
cost of such system, and this last part is again suggested as a continuation to 
this work and not part of it. We will be giving in Chapter 7 our concluding 
remarks. We will summarize the results and discuss future work that we 
suggest necessary before the completion of this project as a whole. 
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Chapter 2 


KINEMATICS AND DYNAMICS OF THE 
CENTRIFUGAL RELAY DESIGN 


Understanding the kinematics of the centrifugal relay system is of great 
importance to this project. Before we start our kinematic study, we will give 
a brief overview of the whole launch system in order for the reader to get 
a clear picture of the concept. We will also discuss design considerations 

for the centrifugal launchers and come up with a conceptual design for the 
rotors. 


2.1 OVERVIEW 

The overall launch system is outlined in Figure - 4 [24]. For the first step of 
the launch process, a large flywheel is geared to a small primary centrifugal 
launcher. The power source could be nuclear, bioconversion, or physical 
conversion of solar energy [19]. Ferromagnetic ball bearings are released near 
the axis onto a curved track, and exit into a guide barrel. The spherical shape 
of the projectiles would minimize the influence of solar winds according to a 
study by Brian Von Herzen [28], although we are not concerned in this thesis 
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Conventional motor-driven 
fly-wheel spun up with a 
solar or nuclear electric system 



V = 0 
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Mechanism 



Relay Payload 
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Figure - 4. Components of a launch system with one centrifugal relay. 
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with the external perturbations from gravitational fields and solar winds. 
For position and timing measurements, each ball interrupts at least two sets 
of horizontal and vertical laser beams directed upon fine-scale photodiode 
arrays. For aiming, sets of top/bottom and left/right coil pairs are energized 
to impart small lateral velocity increments. Solenoidal coils surrounding the 
flight path are used to exactly time the arrival of each ball at a target rotor 
0.1—10 km downstream. 

The stream of balls so launched first encounters a spinning relay rotor. 
In a spaceborne system, deflector shields would protect the relay against 
rare misaimed balls. The rotating balls land on a guide track with velocities 
matched to the track velocity so that they roll up the track. Present plans 
call for the relay rotors to be made of non-ferromagnetic composite materials. 
Coils imbedded near the rotor axis slow down the motion so that, when a ball 
rolls back down the track, it leaves going almost exactly forward or backward. 
Forward-going balls go on to impel the next relay or the payload assembly; 
backward-going balls are collected at the primary launcher. Additional tracks 
on each rotor assembly are used to process balls which land near axis to 
maintain nearly constant rotational energy and/or are being passed backward 
for eventual collection and reuse. 

For a detailed description of centrifugal relays for spacecraft propulsion, 
the reader may refer to “Centrifugal Relays for Spacecraft Propulsion” by 
Clifford Singer and Richard Singer [24]. 

For the purpose of this thesis, we will be considering just one relay and the 
balls are only going backward to the primary launcher. The method is general 
and the results can be generalized easily for forward going balls impelling 
other relays or the payload. Before doing this, however, it is desirable to 
treat the one relay case in all details as completely as possible. We will 
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start here by describing the kinematics involved in this project and show 

how the centrifugal launchers can be used as energy and momentum source 
for spaceborne propulsion systems. 

2.2 KINEMATICS 

The purpose of these kinematics calculations is the determination of momen- 
tum transfer. We will start by studying the case of a one ball launch and see 
how that ball transfers momentum to the payload. We will then generalize 
the results and write a recursive formula that gives the momentum gained by 
the payload as a function of the number of balls launched. This study will 

allow us to estimate the amount of initial stored energy required to launch a 
payload. 

2.2.1 Momentum Transfer after the Launch of the First Ball 

The approach here is similar to the one followed by Clifford Singer and 

Richard Singer [24], The general linear momentum conservation equation 
can be written for this system as: 

mv in + MV in = mv out + MV ouU (j) 

where m is the mass of the ball and M is the mass of the rotor catching the 
ball or. eventually the payload. (?„ is the velocity that the ball has when it 
leaves the primary launcher and tf„, is its velocity after it impacts the relay 
or payload, and rolls back off the track heading back toward the primary 
launcher. K„ is the relay/payload velocity before the launch of this first ball, 
and V ml is its velocity after this first launch is completed. 

The energy conservation equation can be expressed in terms of the rota- 
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tional energy E rot , gained by the relay as: 

E,«= + \MV?\ - + l -MV^ , (2) 

^ V ^ S " ■ S/ ^ 

total initial kinetic energy total final kinetic energy 
or, 

2 E„, = m(vl - + M(Vl - O- (3) 

Let us denote by At; the decrease in ball velocity, 

Av = v in - v out . (4) 

Then, the momentum lost by the ball is mAv. It corresponds to the momen- 
tum gained by the rotor M(V out — K n )- So, the momentum transfer can be 
written as: 

M(VL - V in ) = mAv. (5) 

Equation (3) can be rewritten as: 

2 E rot = m( v in - v out )(v in + v out ) + M(V, n - Kui)(K n + Vout), 
and, using Equations (4, 5), it becomes: 

— #■ -♦ 

2 E rot = mAv(v in + v out ) - mAv(V in + V out ). 

Now, we can still substitute (tf,„ — Av) for v out from Equation (4) and, 

(Vin + jtfAv) for V out from Equation (5) in the equation above and get: 

2 E rot = mAv [2(ui n - V in ) - Au(l + -^r)j . (6) 

In the case of the payload, there is no rotational energy transfer 
(E rot — 0). Equation (6) can then be written as: 

2(v in - Kn) = Au(l + ^)> 
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or 


, using the definition of Av. 


Vj n — Vout i I m 

A 1 \A 


2(t )j n ~~ Vjn) 
m 1 

M 


and finally, 


2^„ + ^«(-l + S) 


i iE 
^ ’ M 


( 7 ) 


We can also determine the increase in payload velocity as: 

V out - V in = ^ An = ^(*« - ^ut). 

Substituting for tw from Equation (7), the increase in payload velocity be- 

comes: „ . m 2(t?j, - Kn) (8) 

Vout-Vin- M ! + ^ 

Now if we go back to the case of a relay, the rotational energy is generally not 
zero’ It can be obtained from Equation (6) provided An has been determine . 
In the next chapter, we will determine the final velocity of the ball, >w a 
we analyze its interaction with the relay using Lagrangian mechanics. e 
have to mention here that our study will assume an infinite mass rotor. In 
practice, the rotor is about 10 kg and the ball is just a few grams. The ratio 
m /M is of order 10“, which justifies our assumption of infinite mass relay 
H we were to take into account the finite mass of the relay, the Lagrangian of 
the system would be more complicated and the angular velocity of the re ay 
would not be constant but would be an extra unknown to be deteruun . 
This would complicate the analysis. We will thus simplify the Lagrangian 
and study the interactions of the ball with an infinite mass rotor, keeping in 
mind that we have an error of the order of (m/M)’ stepwise according to 
Equation (8), where the increase in velocity of the relay would be. 
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Knowing that m/M ~ NT 4 , we believe this is a quite good approximation 
and we will definitely consider in the next chapter an infinite mass rotor. We 
can also estimate at this point the cumulative error after the launch is com- 
pleted. We will be throwing balls on this rotor to increase its velocity until it 
becomes comparable to the velocity of the incoming balls, v tn . According to 
Equation (8) again, the total number of these increments would be M /2m. 
So, the total error can now be determined by multiplying the total number 

of increments by the stepwise error. The cumulative error is estimated to be 
of the order of m/2M. 

In the case of an infinite rotor mass, the rotational energy gained by the 
relay is zero, making the cases of relay and payload similar. Furthermore, if 
we start for this first ball launch with an initial zero velocity for the relay 
(Vin = 0), we get from Equation (7): 


V out ~ “^in, 

and from Equation (8): 

t 7 _ m _ 

K out - 2—v in , 

which represents the limiting case of an elastic collision where the ball trans- 
fers twice its momentum to the relay. Now, we can determine the actual 
rotational energy if we consider a step by step analysis and write the angular 
momentum conservation equation. The angular momentum lost by the ball 
will be gained by the rotor. Knowing v in and v out from the Lagrangian anal- 
ysis, we can write the angular momentum L before and after the interaction 


18 



respectively as: 


Lin = r in X mv in , 

and 

^out r out x 

where r in is the position vector when the ball impacts the track on entry, and 
f out its position at the exit of the track when traveling back to the primary 
launcher. Then, the decrease in the ball angular momentum, L in - L out will 
correspond to the increase in angular momentum of the relay, A L rotor . And 
finally, if we know the moment of inertia of the rotor, I rotor we will determine 
the increase in its angular velocity Aw from: 

— * 

^Lrotor ^rotor . 


2.2.2 Momentum Transfer after the Launch of n Balls 

Let us now consider a sequence of such interactions and derive a recursive for- 
mula giving the increase in relay/payload velocity as a function of the number 
of balls launched. The relay/payload initial velocity, before the launch of the 
first ball, is VJ n . After the impact of the first ball, the velocity becomes: 

l 7 . 1 O m ( V in — K'n) 

* tn T & ~r~z — : — , 

M l+fi 

as shown in Equation (8). We can now write the final velocity of the re- 
lay / payload after the second ball is launched: 


4. O m (Vin Kn) ^' n ^ M ^ lY-St”'* ) 

M 1 + t M l + £ 


+ 2 | 2— jgkzik l _ fo^.1 2 (3n ~ Vin) 


M i + fi 
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We also derived the result for the case of three balls launched (n = 3): 


Vin + 3 


2 m (gn - Kn) 

M 1 + f _ 



We went on and wrote the velocity resulting after 4 balls have been launched, 
and also after 5 balls were launched. We then expressed the general formula 
giving the final velocity of the relay/payload after n balls are launched as: 



We then wrote a little program (cf. APPENDIX), that uses the expres- 
sion derived above and computes the velocity of the relay after n balls are 
launched. We plotted in Figure - 5 the generated data to show how V out 
varies with n. We can see that the velocity of the relay increases almost 
linearly with the number of balls launched, but then, as the relay is going 
faster and faster, the effect of one ball becomes less and less important. Very 
little momentum is transferred when the rotor is moving at speeds similar to 
the speed of the incoming balls, until it reaches the speed of the balls and 
they will not be able to catch up with it. In the practical case considered 
here, some 8000 balls would be needed to reach a velocity v rtlay = 1698 m/s. 
Only the first 10% of the total number of balls needed to reach this limiting 
cas e of a relay moving, practically as fast as the balls, will provide half of 
the final rotor speed while the remaining 90% of the balls will provide the 
remaing half of the total velocity increments. This suggests that one should 
accelerate the relay until it reaches about half the velocity of the incoming 
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Number of balls 


Figure - 5. Velocity of the relay versus the number of balls launched. 
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balls as further acceleration would not be cost effective since the balls are 
not as efficient in terms of momentum transfer. 

2.3 DESIGN CONSIDERATIONS FOR THE CEN- 
TRIFUGAL RELAYS 

One major design consideration is the fact that the relay must process a 
large number of balls quickly and efficiently. The rotor is a disk or a cylinder 
spinning at high angular velocity with guide tracks to process the balls in 
the necessary direction. The balls are either caught near the axis of the rotor 
and launched forward or caught nearer the tip of the rotor, slowed to a stop 
at the axis and then propelled forward or backward using another track. The 
objective of the first type of processing is to increase the velocity of the balls. 
The second type of processing increases the angular and linear momentum 
of the relay. For the determination of the exit velocity and the forces of 
constraint on the ball, a track shape needs to be chosen. If the ball starts at 
the center of the rotor and is launched tangentially, it leaves with twice the 
tip speed of the rotor, where the maximum tip speed depends on material 
properties of the rotor [23]. The shape analyzed here is the logarithmic spiral. 

The objective is to develop a method that can be used to analyze different 
shapes quickly and easily. 

Original studies of the kinematics of ball/rotor interaction idealized the 
relatively small balls as point masses [24]. Then an equation of motion treat- 
ing the pellet as a ball of finite radius was derived [29]. One constraint to 
be satisfied is that the balls must not bounce nor slide during the catching 
process to insure a maximum exchange of momentum between the ball and 
the rotor. This constraint is satisfied if the velocity vectors of the point of 
impact on the ball and the blade are equal. We have started detailed kine- 
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matic studies by assuming that the velocity vector at the edge of a ball rolled 
off a previous rotor must match the velocity vector of the contact point on 
the next rotor it encounters. This condition is met by slightly adjusting the 
aiming and timing so each ball launched encounters the next rotor at the 
correct distance from this rotor’s axis. This removes a degree of freedom 
which complicated the analysis made in the paper “Centrifugal Relays for 
Spacecraft Propulsion” [24], This should considerably simplify the next step 
in the kinematic analysis-writing a computer program which follows the tra- 
jectories of all objects in the system in the presence of a gravitational field. 

To obtain an even more complete description of ball/rotor interactions, 
AUTOLEV program may be used. This program, by Schaechter and Levin- 
son [30], automatically generates FORTRAN coding for complete solution for 
rigid body kinematics of finite mass balls interacting with a free-flying rotor 
in a gravitational field. While including the effects of finite ball/rotor mass 
ratio is not expected to be important for the kinematics of the overall mission 
performance, it will determine the exact amount of magnetic drag needed to 
hold the balls near the rotor axis long enough that they will roll back off 
at the right time. The same method can also be used to examine various 
schemes for emplacing and activating these coils. To date, the AUTOLEV 
program has been successfully used for balls landing with properly matched 
velocities, rolling up towards the axis of a rotating relay, and rolling back off 
the track for relaunch [31]. The program is designed so that investigation 
of additional effects like small velocity mismatches at catch, finite ball/rotor 
mass ratios, and gravitational fields are straightforwardly accomodated by 
changes in the program input. 
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2.4 CONCEPTUAL DESIGN OF THE RELAYS 

The first task in designing the centrifugal relays is coming up with a con- 
ceptual idea. The rotor has to have guide tracks for the peUets. It should 
spin at very high angular velocity processing balls in both directions. The 
rotor should include guidance tubes in the path of the incoming pellets for 
position control. A protection shield has to be placed in front of the guide 
tubes and rotor to prevent errant balls from impacting these systems. A 
supporting truss structure is needed to mate the rotor, the guidance tubes 
and the protection shield. To minimize the total system mass, the size and 
the mass of the whole relay may need to be kept small. This conceptual 
design can be summarized schematically as illustrated in Figure - 6 [32]. 

2.5 DESIGN CONSIDERATIONS FOR THE RO- 
TOR 

The basic concept of a centrifugal projectile accelerator is that a projectile 
can be accelerated centrifugally by either being constained to move along a 
channel in a rotating disk or arbor, or by being released from the edge [23], 
Figure - 7 illustrates such launcher [23], 

The rotor is a disk or a cylinder spinning at high angular velocity with 
guide tracks to process the balls in the necessary direction. The “forward” 
or “backward” launching depends on the objective of the launch. It could be 
increasing the velocity of the ball, or the angular and linear momentum of 
the relay. The idea in this design is to have one rotor accomplish one task. 
For example one rotor would catch the balls near the tip and slow them to 
a stop at its axis increasing the rotational and transitional velocity of the 
relay. The balls would be collected at the axis of this rotor and then fed to 
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Figure - 6. Elevation and plan view of a centrifugal relay. 
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another disk to be propelled in the necessary direction (forward or backward). 
A display of this concept is given in Figure - 8 [32]. The disks are all rotating 
at the same rate. The balls are fed through the shafts to the correct track. 

In practice, instead of having several disks, we will just use a multiple track 
rotor. 

Another point discussed here is whether the rotor should be a cylinder 
with a track bored through it or a rotating disk with a surface track. The 
latter case may seem more interesting for reasons of machining and catching 
incoming balls. It is easier to catch the ball on a surface track then in a small 
hole. But the track bored through a cylinder may be more practical for the 
transfer of the ball from one track to the other at the center. Conceptually, 
the best design would be a combination of these two alternatives. We would 
have a surface track starting at the edge of the rotor and continuing into a 
path bored through the rotor. The track will end at the axis of the cylinder. 
The second track will start from the center at a lower level and emmerge to 
the surface toward the tip of the rotor. 

Typical materials of interest for use in a pellet launcher are given in 
Table - 1 [33], where V c represents the critical velocity, p is the density and 
a the yield strength of the rotor material. To build a launcher that ejects 
pellets with velocities less than V c , we can just use a tube of uniform cross 
section. Above V c , the pellet launcher must be tapered for part of its length, 
increasing in thickness toward the axis [33]. 

The tube may be subject to considerable wear due to friction and abrasion 
from the pellets, and should be designed for easy maintenance. This could be 
accomplished by providing the tube with a removable liner, and by designing 
other high - wear parts for easy removal and replacement [33], 

The U. S. Department of Energy has designated Lawrence Livermore 
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Figure - 8. Multiple track relay concept. 
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Table - 1. Materials of interest for 
a rotary pellet launcher 
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National Laboratory the lead laboratory in its flywheel rotor and technology 
development program. Although the objective of that program was storing 
energy, the rotors that were developed can be modified to launch projec- 
tiles [23]. One of the near-term objectives of this program is to develop an 
efficient, economical, and practical composite flywheel. State-of-the-art com- 
posite rotors developed in the program [34] have maximum tip speeds < 1 
km/s. Their radii are typically 10 — 20 cm, and masses are a few kilograms. 
These centrifugal launchers are summarized in Table - 2 [23], 

For the determination of the exit velocity and the forces of constraint 
on the ball, a track shape needs to be chosen. The launch should be such 
that the exit velocity is approximately twice the tip speed of the rotor. The 
maximum tip speed, V c is proportional to (cr/p) 1 ^ 2 [23]. 

The shape considered is the logarithmic spiral. An illustration of the 
track shape is given in Figure - 9. The method developed should be general 
and could be applied to analyze different shapes quickly and easily. This is 
the subject of the next chapter, a detailed analysis of the launching process. 
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Rotor 

State-of- 

the-Art 

Near- Term- DOE 
Objective 

Demo 

Maximum energy 
density (MJ/kg) 

0.24 

0.32 

0.57 

Maximum rotational 
energy (MJ) 

0.62 

3.6 

3.4 

Mass (kg) 

2.6 

11.2 

6 

Radius (cm) 

20.3 

30.5 

120 

Maximum frequency 
(rpm) 

49,000 

35,400 

12,000 

Maximum tip 
speed (km/s) 

1.04 

1.13 

1.5 


Projectile 

State-of- 

the-Art 

Near-Term-DOE 

Objective 

Demo 

Mass (g) 

26 

112 

60 

Maximum velocity 
(km/s) 

1.04 

1.13 

3 

Maximum energy 
(MJ) 

0.014 

0.072 

0.27 


Table - 2. Rotors and launching capabilities 
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Chapter 3 


ANALYSIS OF THE LAUNCHING 

PROCESS 


Lagrange equations with multipliers are used to analyze the launch pro- 
cess. The track shape is represented as a constraint equation between the 
radial and angular positions. A second constraint equation represents the 
fact that the ball is restricted to roll on the track. We wish to avoid skidding 
because frictional losses would make the system non-conservative. 

The velocity of the center of mass of the ball along the track is: 


r -t- r'tr, 




and 


Ht) — + <*>t + C(0> ( 11 ) 

where r is the radial displacement of the ball, u> the rotation speed of the 
rotor, C the angle that defines the position of the ball along the track, and 

*° any initiaI an ® le that might be present as illustrated in Figure - 10. If we 
let m be the mass of the ball, then the kinetic energy is: 

T = \ m + r > + C) 2 ] + (12) 
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where J is the moment of inertia of the ball defined as 


J =5^’ 


( 13 ) 


1 being the radius of the bal1 and « the angular velocity of the ball around 
its center of mass. These quantities are all defined in Figure - 11. 

The Lagrangian is: 

L = T-V = T (14) 

since there is no potential energy. 

The constraints equations are: 


fi = r ° e “ at) - r{t) = 0 , ( 15 ) 

and 


f 2 = £ a (t) - s(t) = 0, (16) 

where s, the track shape length defined in Figure - 11, is related to r as 
follows: 


ds 2 = dr 2 + r 2 d( 2 , 


°r> 



(17) 

(18) 


Now, from Equation (15), we can write: 


Then, 


C = — In — 

a r 0 


dr 


1 


ar 



= 1 + -T, 


(19) 

( 20 ) 

( 21 ) 
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Figure -11. Rolling of the ball on the track* 


or, 



( 22 ) 


Integrating this equation and using the initial condition 5 = 0 when r = r 0 
to determine the constant of integration, we can write s(t) as a function of 
r(t): 

(/h^) (r -r°), (23) 

and / 2 becomes: 

/ 2 = ta(t) - (r(i) - r 0 ) = 0. (24) 

The equations of motion are: 

i (<>£\ df* 

dt \ dr J dr 1 dr 2 dr ’ 

± W _ d ± - X d/i ^ . dh 

dt \dc) d( 1 d( + 2 dC 

and 

± (<Hl\ _ ^ _ i d f* . i 

dt \ dd y da 1 da 2 da 

Dividing Equation (25) by m gives: 

f ~ r(u + C) 2 + — + — Jl + ^ = °- (28) 

m m V a 2 

Equation (26) can be written as: 

r 2 ( + 2rr(u> + () - ar— = 0. (29) 


(25) 

(26) 

(27) 
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~ *2, 

j' =0 - 


(30) 


Dividing Equation (27) by T gi 


From Equation (30) combined with Equations (16, 23) differentiated twice 
with respect to time, we get: 



Solving for A! from Equation (29) we get: 


Aj 2 lo , 2 r - 

m = T r + a ^ + a C ' ( 32 ) 

Now, substituting for Aj from above and A 2 from Equation (31), and using 
the following set of equations derived from Equation (10) 


and 


Equation (28) becomes: 


r i_ _2_ 

L a 2 r a 2 r 


or simply, 




1 

a 2 r 


r fl + 

1 ( 

1 \ 

1 I 

ml 2 l 1 + a 2 ) 

+ ?] + 

JH 

■ 2u> 

2a; 1 


1 — 

a a J 

ru 2 = 0, 


f - K 2 r = 0, 


where 

, .2 

2 _ W 

(' + i) (i + <?)’ 


(33) 

(34) 


(35) 

(36) 

(37) 
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and 


3 = 


We can also write K 2 as: 


I < 2 = 


mP 


5 u> 2 


(38) 


7 (* + ?) ’ 

taking into account the fact that the projectiles have a spherical shape, which 
means: 

*-!• 

The corresponding solution to the governing equation is: 

r(<) = A cosh K(t - t 0 ) + 5sinhA'(< - t 0 ). 

Differentiating this equation, we can write: 

r(<) = AK sinh K(t - t 0 ) + BK coshA^f - f 0 )- 
Now, at t = < 0 j when the ball is released near the center, we have: 

r(t o) = r 0 , 

and 

r(<o) = 0. 

Consequently, 


and 


A = r 0 , 


B = 0. 


The trajectory of the ball on the track is then given by: 

r (0 = r 0 cosh K(t — t 0 ). 


(39) 


(40) 


(41) 
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One quantity of interest is the exit velocity. It can now be determined, 
and we can study its behavior as a function of the design parameters such 
as the track shape constant a. 


3.1 EXIT VELOCITY 

The position of the ball on the rotating track with respect to a fixed frame, 
with 0 O = 0, is given at all times by: 

^)=f roe “ < "' COSM+C( ‘ )) 'l (42) 

{ sin(u,t + «t)) J ' 


We can get the velocity vector by differentiating the position vector above, 


( r 0 oCe a<:(t) cos (wf + £( t )) - r 0 (u> + C)e a «‘> sin(u;< + £(*)) 
\ r 0 aCe a <(‘) sin (art + £(<)) + r 0 (w + C)e a «‘) cos (ut + ((t)) 


(43) 


Then, 

v 2 = rle 2a< [(£ + u ) 2 + a 2 C 2 j 
Using Equation (33), v 2 can be written as: 


„2 __ 2 
v = r 


( r N 2 
— + u 
\ar 




(44) 


(45) 


Now, from Equation (41) we get: 


r = Kr 0 sinh K(t — f 0 ). (46) 

This can be combined with Equation (41) to express a relationship between 
r and r, 

r = Kyjr*-rl. (47) 
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Substituting for r from Equation (47) and K from Equation (38) in Equa- 
tion (45) we get: 


v 2 = rV 


And finally, we can write the exit velocity v exit as: 


1 +a 2 


(48) 



where v tip is equal to ur x . 

Also from the geometry, as outlined on Figure - 12, we can write: 


tan 7 = ^ = M^ 

dx dx/d(' 

where tan 7 is the slope of the tangent to the track at the exit, and x and y 
are given by: 

x = e a< cos(( + ut), 

and 

y = e° c sin(C + ut). 

We can also write tan 7 as: 

tan 7 = - Sm ^ - ~ l ~ cos (C + 

a cos(C -f- ut) — sin(£ + ut) 


Now, let us express C in terms of 7. We can write, after inspecting Figure - 12: 


Tangent to circle 2 ^ + ^ ^ (rc - ^)] 



Figure - 12. Angular information of the ball at the exit 
from the primary launcher. 
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where 0 is the angle between the tangent to the rotor rim and the tangent 

to the track at the exit. Using this relationship between the angles, tan 7 
becomes: 


fggCfl + 7) + sin(/? + 7) _ sin 7 
a sin(/? + 7) + cos(/? + 7) — cos 7 

This equation can be rewritten as: 


a cos (0 + 7) cos 7 + sin(/? + 7) cos 7 = a sin (0 + 7) sin 7 + cos (0 + 7) sin 7, 


or simply, 


-a cos /3 = - sin/?. 


Thus, the main result of this geometrical derivation is: 


tan 0 = a 


or, 


cos 0 = 


1 

VTT^' 


If the ball leaves from the center (r 0 = 0) tangentially (/? = 0) Equa- 
tion (49) becomes: 

Vexit = 

and we can see that the ball leaves with almost twice the tip speed of the 
launcher. Felber [23] stated that the ball leaves with twice the tip speed in 
the same conditions of launch. But in his case, the ball was considered as a 
point mass with zero inertia (t -+ 0), which means that all the energy gained 
by the ball is used for the translational displacement of the ball, whereas 
in our case, some of the energy is used to spin the ball around its center 
of mass. So, the resulting translational velocity is less than in the case of 
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Felber’s calculations, and in general, the exit velocity does depend on the 
mertia of the projectile. We can actually find his result from our general 
expression just by setting J = 0 (no inertia) in Equation (37) which gives: 




ur 


i + £' 


Then, substituting for r from Equation (47) and K from above in Equa 
tion (45), we find: 


Vcxit — Vtip 

2-(- 


-Y 1 1 



*1' V Vf 

'J VTT^ 


or, as a function of the angle /? at launch, 


Vextt — Vu p 




which is what Felber [23] published. In this particular case, we can see that 

when the ball leaves from the center (r„ _ 0) tangentially (0 = 0), it leaves 
with twice the tip speed (u e;ti< = 2v tip ). 

The dependency of the exit velocity on the track shape constant a is 
illustrated in Figure - 13. 

Another important point about this launch process analysis is the deter- 
mination of the time delay needed to insure that the ball is launched exacly 
forward. The ball will be kept near the rotor axis while launcher is rotating 
until the track is oriented in such a way that when the ball rolls off the edge 
it is going exactly forward. 
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a (track shape constant) 


Figure - 13. Exit velocity of the ball versus the track 

shape constant. 


45 



3.2 TIME DELAY NEEDED TO INSURE LAUNCH 
EXACTLY FORWARD 

The position of the ball on the rotating track is given by Equation (42) and 
its velocity by Equation (43). In order for the ball to go exactly forward 
(horizontally), v at the exit should be of the form: 



Let t x be the time at which the ball leaves the launcher, and t 0 the time 

at which the ball is released near the center of the rotor. The origin of 

time ( t = 0), is taken to be at the beginning of the spinning process of the 

launcher. We can determine the time it takes for the ball to roll down the 

track all the way to the edge, which is (t 1 - t 0 ), by solving Equation (41) at 

the exit (t = t x ), where r{t ) will be set equal to r,. Then we can get the total 

launch time t x by setting Equation (43) at the exit equal to Equation (50). 

Finally, the delay time t 0 can be determined just by combining these two 
results. 


3.2.1 Determination of (< x - * 0 ) 

At the exit, Equation (41) becomes: 

r (*i) = n = r 0 cosh K(t x - t 0 ). ( 5 i) 

This equation can be solved for (t x - t 0 ) as: 

<'-<0 = j?cosh- (^)- (52) 
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3.2.2 Determination of f x 

Equation (43) can be rewritten as: 

v = | cos ( a '^ + 0 ~ r(u> + C) sin(u;t + Q 
\ sin(wt + C) + r(u> + C) cos(u>t + C) 


v H 

Vv 



( 53 ) 


At the exit, t — f x , the two components of the velocity vector have to 
satisfy the condition expressed in Equation (50). We must have w = 0 or, 

rj sin(u>*i + Ci) + r x (uj -f C x ) cos(o;t x + C x ) = 0, (54) 

where the subscript u l” means that the quantity is taken at the exit. This 
equation can also be written as: 


r i tan(wij + Ci) + r x u> -f- r x — — = 0. 

ar i 

Then, 

+ Ci = arctan f r-w] + kir, 

L a r x J 


(55) 

(56) 


where r x can be obtained from Equation (47) and Ci from Equation (19) 
where we replace r by r x . Finally, we can write the expression of t x as a 
function of the parameters of the system: 


t\ 


u 


^arctan 




(57) 


The parameter k is adjusted so that the results are physically acceptable. 
The value of t x must be greater than (i x - t 0 ) found previously. We can 
check the accuracy of the results here by plugging the expression of < x above 
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into the first component of the velocity vector at the exit, vjf and it should 
satisfy the equation Vfj = v erit . 

These calculations can now be applied to the practical case we are con- 
sidering and the actual values of t 0 and t\ can be determined. This will allow 
us to have a good picture of the launching process. 

3.3 APPLICATION 

The practical case considered is as follows: 

• The radius of the rotor is: r x = 0.1 m. 

• The distance from the center of the rotor at which the ball starts the 
acceleration process is: r 0 = 0.01 m. 

• The track shape constant is taken as a = 1. 

• The launcher is rotating at an angular velocity u = 10000 rad/s, 
which gives a tip speed of u ftp = 1 km/s and a revolution period of 
T = 2 7r / u = 0.63 ms. 

The results are as follows: 

1. The ball leaves the launcher with an exit velocity v ex n = 1.7 km/s. 

2. The timing is such that: 

• *i — = 0.5 ms from Equation (52), 

• = 0.9 ms from Equation (57), where k = 4, 

• and t 0 = 0.4 ms. 
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The launch process for one ball lasts 0.9 ms, which is about one and a 
half revolutions of the rotor. For the first 0.4 ms of the launch process, the 
ball is kept near the axis of the rotor at the distance r 0 , then it is released. 
After another 0.5 ms, it will roll off the track horizontally at a velocity of 1.7 
km/s. We also substituted for t\ in t;// at the exit using the numerical values 
above to check their consistency and we found v H = v exit = 1.7 km/s, which 
confirms the accuracy of the results. 

The launch process is illustrated in Figure - 14 through Figure - 16 where 
we display the position of the track at different sequences of the process. 

The last part of this launch process study consists of the roll versus skid 
analysis of the ball on the track. This will be the subject of the final section 
of this chapter. 

3.4 ROLL VERSUS SKID ANALYSIS 

The fundamental assumption for this system is that it is conservative. How- 
ever, the high speed at which the ball strikes the track on entry and the 
acceleration at which it is flung outwards may cause the ball to skid on 
the track. Skidding would cause frictional losses making the system non- 
conservative. S. Pang considered this and drew some conclusions about the 
sliding tendency and the system parameters [29]. Skidding is considered here 
in the context of a high vacuum environment where the coefficients of friction 
are largest, as illustrated in Table - 3 [29]. A higher coefficient of friction 
between the ball and the track provides a higher tangential force, raising the 
critical point at which the ball begins to skid. However, if skidding indeed 
occurs, the frictional losses are severe. 

The purpose of this analysis is to investigate the conditions that cause 
skidding and to determine design parameters to avoid skidding. In addition 
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Figure - 15. Track orientation at the instant t = t 0 = 0.4 ms. 

The ball is now released near the center of the lanucher. 
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Figure - 16 . Track orientation at the instant t = <, 
The ball now rolls of the track. 


J 


0.9 ms. 
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Material Pair 

In Air 

In High Vacuum 

Al-Al 

0.78 

1.57 

Be-Cu-Be-Cu 

0.58 

• 

Brass-Brass 

0.32 

1 

Copper-Copper 

1.04 

I 

Be-Cu-Brass 

0.38 

0.9 

Cu-Steel 

0.55 

0.96 

Cr-Cr 

0.85 

1.30 

Stainless Steel- 

0.51 

0.93 


-Stailess Steel 


Table - 3. Table of coefficients of friction 
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to the normal force Xi restricting the ball to move in a spiral, a tangential 
force A 2 constraining the ball to roll on the track is taken into account, as 
explained in the beginning of this chapter. The frictional force available 
between the two surfaces is y\ u y being the coefficient of friction. If this 
available friction is less than the required friction A 2 , skidding will occur. So, 
to check for skidding, we need to check if A 2 > y\ x or if the ratio X 2 /X 1 is 
greater than y. The methodology followed here is as described in [29]. To 
determine A 2 and A! we just go back to the beginning of this chapter and 
solve Lagrange’s equations for these two forces. 

The results are: 



from Equation (31). So, the tangential force per unit mass acting on the ball 
can be written as: 


X 2 _ 2 ( au 2 \ 

m ~ 1 yjT+^) roCOShK{t ~ to) ' (»») 

using Equations (13, 38) and the expression of r(<) from Equation (41). 
Now, from Equation (32) combined with Equations (33, 34), we can write 
the expression of the normal force per unit mass acting on the ball as: 

A^ _ r r 2 2 u> . 

m a 2 + a 2 r + ~a~ r ' 

Again from the expression of r(t) in Equation (41), Aj/m becomes: 
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Ai 5 r 0 u> 2 

m ~ 7 1 + a 2 l cos h -^(t - t 0 ) + sinhA'(i - < 0 )tanh A"(f - t 0 )+ 


•&< ~ a 2 sinh/f(/ - t 0 )j 


( 59 ) 


Now, we can write the ratio of the tangential force to the normal force 


as: 


2a\/ 1 + a 2 

Al 5 [1 + tanh 2 K(t - < 0 )J + 2^35 v/Tf ^ 2 tanh K(t - (60) 


Several calculations were done to see how the design parameters affect 
the sliding tendency of the ball. We can see that the mass and the radius 
of the ball have no effect on the sliding tendency. We can also see that, as 
the ball leaves the center (when t = t 0 ), the ratio is independent of K, since 
tanh K(t - to) = 0. Therefore, the sliding tendency due to the choice of 
u> (or K), is maximum at the end of the track and not at the beginning. 
At that instant, t - t 0 = 0.5 ms, as determined previously in the practical 
case considered. So, we plotted A 2 /A, versus w at the exit of the track. 
This is illustrated in Figure - 17 and we can see that the ratio is practically 
constant. This can be understood from the fact that the Hyperbolic tangent 
varies between 0 and 1, and for the values of u> that we are interested in 
(around 10,000 rad/s), varying w does not make the hyperbolic tangent vary 
much. The denominator in Equation (60) is then almost contant making 
the ratio independent of K, which, in turns is proportional to the angular 
velocity «. This implies the irrelevency of w as a design parameter in this 

analysis. Figure - 17 illustrates the ratio A 2 /Ai versus u for 4 different values 
of the 
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track shape constant a, and we can see that the sliding tendency increases 
with increasing a. So, the design parameter affecting the sliding tendency 
is the track shape. We then maximize the ratio as a function of time to see 
how the geometry affects the sliding tendency in the worst case. We know 
that the argument K(t - t 0 ) is positive. This makes tanh K{t - t 0 ) vary 
between 0 an d 1, an d then, the denominator in Equation (60) between 5 and 
(10 + 2 v / 35\/l + a 2 ). Consequently, the ratio is such that: 

- - 2 qv / 1 + q2 . A 2 „ 2aVTT^ (\ 2 \ 

V A i/ m ,n io + 2vwrr^- a, - 5 ~ l XTy 

' 1 / max 

And finally, we give in Figure - 18 an illustration of the relationship between 
the maximized sliding tendency (A 2 /Ax) ma;r and the track shape constant a. 
Flatter track produces a greater tendency to skidding. Once the materials 
are chosen, p is fixed and the track shape constant a can be chosen not to 
exeed a critical value. It has to be such that A 2 /A a ( required coefficient of 

friction) is not greater than ft (available coefficient of friction). Otherwise, 
skidding will occur. 
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Figiire - 18. Required coefficient of friction versus 
the track shape constant. 
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Chapter 4 


ANALYSIS OF THE CATCHING 
PROCESS 


The ball, launched from the primary centrifuge according to the process 
described in detail in the previous chapter, is then caught on the second 
centrifuge. The resulting momentum transfer is used to propel this centrifuge. 

4.1 CATCHING CONDITION 

For a maximum momentum exchange, the velocity vectors of the points of 
impact on the ball and the rotating blade must be equal to insure no bounc- 
ing, nor sliding of the ball on the track. The normal components of the 
two velocity vectors have to match so that the ball does not bounce off the 
track, and the tangential components have to be equal to insure no sliding 
of the ball on the track. This catching condition can be summarized by the 
following equation: 


V\ — = 


(61) 
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where, tq is the velocity of the tip of the ball, and v 2 the velocity of the 
impact point on the track where the ball is to land. Solving Equation (61) 
will allow us to determine the exact position and timing for catching. 

4.2 DETERMINATION OF CATCHING POSITION 
AND TIME 

At the exit of the primary launcher, the center of mass of the ball has a linear 
velocity v ei jt, as mentioned in the previous chapter. The ball is also rotating 
at a rate or as explained before. This makes the impact point on the ball 
have a velocity u"! such that: 


vi = 


Vexit - U COS £ 
— u sin £ 


(62) 


where the angle £ is defined in Figure - 18 and u is the rotating speed, 



from Equation (24) in the previous chapter, differentiated once with respect 
to time. This can be, equivalently written as: 



Substituting for r x from Equation (47), also in the previous chapter, taken 
at the exit of the track, we can write: 
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The velocity of the point of impact on the rotating track is: 

~ _ f v re i ay — r 0 u;e af sin(u;< + £) 
y r 0 we° * cos (ut + £) 

where u re/ay is the linear velocity at which the relay is moving. This velocity 
vector is obtained by differentiating with respect to time the position vector 
of the impact point on the blade, 



( v reiay t + roe 0 *" cos(o;t + £) 
roe at * sin(u>tf + () 


( 66 ) 


We have to keep in mind while differentiating that we are only interested in 
one point, the impact point. This means that ( is not a variable here. It has 
a specific value that we have to solve for 
From Equations (61, 62, 64) we get: 

Vrelative ~ U COS £ = —UJroe^' sin (u/t,- + £), 

and 

- usin^ = u>r 0 e“ Ct cos (ut { + £), 
where v relative is the relative speed of the two moving bodies (ball and relay), 

Vrelative = V exit - V relay , (69) 

and ti is the time at which the ball impacts the track. 

When we solve Equations (67, 68) the catching parameters, position on 
the blade, £ and time at which catching occurs, will be determined. 


(67) 

( 68 ) 
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From Figure - 19, we can see that: 


tanf = *L = 

dx dx/d(' 


where 


and 


Then, 


x — r^e ^ cos(C + wf), 

y = r 0 e aC sin(C + tut). 
t.-t- l+atan(C+wt) 

cl tan(£ -(- c ot} 


Let us now consider the following identities: 


(70) 


and 


cos £ = ± 

yl + tan 2 £ 


sin^ = ±— 

y 1 + cot 2 £ 


Using the relationship between (tan£) and (tan ((+u>t)) expressed in Equa- 
tion (70), (cos £) and (sin£) can be rewritten as: 


cos£ = ± — t ~ tan ^ + h?< )) 

v/(l + a 2 )(l + tan 2 (£ -f lj<)) ' 


and 
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Point of impact 
on the blade 



Figure - 19. Catching relay and track. 
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sin£ = j (l + qtan(C + u;0) 

\/(l + a 2 )(l + tan 2 (C + u;t)) 

Now, substituting for (cos£) and (sin£) as above in Equations (67, 68), we 
get: 


n(q - tan(C,- + a;*,)) 

v /(l + a 2 )(l + tan 2 (C t +^)) 


^relative T ~ = v, ‘ == = —r 0 u>e a< ‘ i sin(u>t, + £,), (71) 


and 


u(l + a tan(£,- + wf,)) . 

T ~f = r Q <jje a t' cos(u + £,■). 

\/(l + ° 2 )(1 + tan 2 (Ci + o>t,)) 


(72) 


Dividing Equation (71) by Equation (72) we get: 


T ^r e /at,t, eV /(l + q 2 )(1 + tan 2 (C,- + ^ j) a - tan(£, + wt,-) 
u(l +atan(C,- + u><,)) + 1+ a tan(£, + uti) 


= -tan(C, + wf,). (73) 

This can also be written as: 

-p ^relative ~b q2 )(1 d~ tan 2 ((~, -j- u>f,-)) a(l + tan 2 (£,- -f- ujti )) 

u 1 + a tan(£, + wf,) 1 + a tan(£,- + art,) ’ 

or simply, 


=F 


^ relative a/ 1 d* O 2 

U a 


\/l + tan 2 (£, + u>f,). 


(74) 
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If we square this equation, we can solve for tan(£, +u;<,) as 


tan 2 (C, +<*,*,) = (^l^) 2 (l + 1) _ L 


( 75 ) 


Now, we can go back to Equation (72) and rewrite it using the identity: 


cos(C,- + uti) = + 

v 1 + tan 2 (C, + cot,) 


The equation then becomes: 


u(l -)- o tan(^,' -)- ujf, - )) 

VI + a 2 


= rowe 0 ^, 


(76) 


or, if we want to solve for 


C« = - In 

a 


u(l + atan(£ + u^ ,)) 

r 0 u;Vl + a 2 


1- 


(77) 


where, tan(C, + “>*,) is given by Equation (75). The solution, gives the 
position of the ball on the track at catching. We can also solve for the 
argument (ujt { + C.) of the tangent from Equation (75), and finally deduce t„ 
the time at which catching must occur. This time can be written as: 



UJ 


arctan 



- (i + ** , 


(78) 


where, k is chosen so that the catching time is greater that the time the ball 
exits the launcher. 
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4.3 NUMERICAL APPLICATION 

Let us now consider again the practical system proposed in the previous 
chapter and determine the position and time of catching. Let us assume at 
this point that the second centrifuge is not moving yet. It is just rotating at 
the same angular rate u> as the primary launcher. Consequently, 


^relative = V exit = 1.7 km/s. 


We can also have from the previous chapter (Equation (47)) the value of r 
at the exit, which is, rj = 595 m/s for the particular example analyzed in 
the previous chapter. Substituting for r, in Equation (63), we can determine 
the value of u as: 


We also have: 


u = 841 m/s. 


• r 0 = 0.01 m, 

• and u = 10000 rad/s, as before. 

• The track shape is the same as for the launcher case but its orientation 
is different, and so the constant has a different value. For the catcher, 
the track shape constant is a = -1. 

After substituting for u, v re i ative , r 0 , u>, and a in Equation (75) we get: 

tan 2 (Ci + uti) = 7.17, 


or, 


tan(0 + u/ti) = ±2.68. 
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Now, substituting this result in Equation (77), we fing the catching position: 

C» = —2.3 rad, 

which represents the edge of the track. This result confirms our intuitive 
feeling that the ball had to be caught on the edge of the track to satisfy 

the velocities matching as explained before. Finally, we can solve for U from 
Equation (78): 

ti = 1.3 ms, 

where k was taken equal to 3. 

Knowing the total time of the launch sequence = 0.9 ms from the 
previous chapter, we can deduce the time of flight of the ball between the 
two centrifuges. This time is: 

t flight ~ ti — ti = 0.4 ms. 

And knowing the velocity at which the ball is travelling, we can deter- 
mine what distance, d , initially we should place the relay from the primary 
launcher: 

^ = v exit t flight = 0.68 m. 

etter to start with an initial separation of about 500 m to insure that 
the time the ball spends on the relay is very small compared to the total time 
of flight of the ball for kinematics studies purposes as we shall see later. The 
ball will then be caught several revolution later than the time f, calculated 
above. We have to choose k such that the time of flight will correspond to the 
desired distance. This fact will not be taken into account in this chapter. We 
will consider the minimal separation between the launcher and the relay of 

0.68 m, but we will make the necessary changes in the generalization chapter 
for multiple balls launching. 
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To conclude this section on the catching conditions, we illustrate the 
orientation of the track both at t = 0 and at t = t { in Figure - 20 and 
Figure - 21. 

4.4 BEHAVIOR OF THE BALL ON THE CATCHER 

Now that the first ball is caught smoothly at the edge of the blade without 
any bouncing or sliding, we have to determine its trajectory on the track. 
As for the launch process, Lagrange equations with multipliers are used to 
analyze the behavior of the ball on the centrifuge. The same constraints as 
before are on the system. The first specifies the track shape and the second 
insures that the system is conservative by constraining the ball to roll on the 
track. Equations (10) through (21) from Chapter 3 still apply, although the 
orientation of the angle < and the track shape constant a are different in this 

case. The reader should refer to Figure - 19 where an illustration of the track 
on the catching rotor was given. 

Equation (21) from the previous chapter gives in this case: 

d3 = ~ (j 1 + dr ■ (79) 

This equation can be integrated and, along with the initial condition, s = 0 

when r = r u it gives a relationship between the track length and the radial 
position of the ball: 

s W = ( y 1 + ^2 j ( r i -r). (80) 

Then, the second constraint equation becomes: 

/ 2 = £a(t) + + — j ( r (t) - ri ) = q. (81) 
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Figure - 20. Orientation of the track on the relay at the 

instant t = 0. 
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igure - 21. Orientation of the track on the relay at the 
instant t = t { = 1.3 ms. 

The ball now impacts the relay. 
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The equations of motion turned out to be exactly the same as for the launch 
process, and the governing equation is: 


r(t) = A cosh K(t - U) + £sinhiif(f - *,-). ( 82 ) 

The initial conditions at catching are known in terms of position and velocity: 

r (<«) = r u 

and 


V (^«) — V cxity 

for the center of mass. 

Now, from »(<<) we can get r(t t ). We have: 


v 2 = r 2 + r 2 ( u + 


from Equations (10, 11) of the previous chapter. We can also use the rela- 
tionship between f and r as expressed in Equation (33) in Chapter S and 
the previous equation becomes: 


or. 


112 (' + ?) + \ T ™ + 


( 83 ) 


r 2 ui 2 - v 2 = 0. 


The solution to this equation is: 


• -lru±SD 


r ~ 


where D is defined as: 


f 1 + £) 


( 84 ) 




( 85 ) 
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Equation (84) can be expressed as: 


—aru> ± -(-a 2 ) — a 2 r 2 w 2 

1 + a 2 


At t = f, we get: 


HU) 


— aria> + -f a 2 ) — q 2 r?u> 2 

1 +a 2 


( 86 ) 


(87) 


where the “+” sign was chosen over the sign because r(t t ) has to be 
negative since r(t ) is decreasing from r X) radius of the rotor, to r 0 near the 
center. We also know that the track shape constant a is, in the case of the 
catcher, negative. It turns out that: 


HU) = ~r u 

where rj is the value of r at the end of the launch sequence (at time t x ). 
Now, from Equation (82) we can write: 


r(U) = A. 


Differentiating Equation (82) with respect to time and setting t = t { we can 
also write: 

HU) = BK. 

Identifying these two equations with the initial conditions at catching as 
expressed before we get: 


A = r u 
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and 



The trajectory of the ball on the catcher blade is given at all times by: 


r{t) = r x cosh K(t - t,) - £ sinh K(t - t { ). (88) 

Let us denote the time it takes for the ball to reach the center of the rotor 
by t 2 . Solving Equation (88) at the center ( r(t ) = r 0 ) gives the value of the 
interval of time (t 2 — <,). This interval of time is the time actually spent by 
the ball on the track of the catching rotor while t 2 is the total time from the 
very beginning of the launching process ( t = 0). 

The ball is now going to be transferred to another track on the catcher 
and then be propelled forward or backward. For the purpose of this partial 
study, we will only consider throwing it back to the primary launcher. 

4.5 RETURN OF THE BALL TO THE PRIMARY 
LAUNCHER 

The ball is now on the second track of the catcher to be propelled backward. 
The track shape is always the same and the orientation is opposite to the 
track where the ball is caught making it similar to the track where the ball 
is accelerated on the primary launcher. The trajectory of the ball is then 
known. We just need to replace t 0 in Equation (41) ( Chapter 3 ) by t 3 and 
we can write: 

r (0 = r 0 cosh/f(< - * 3 ), (89) 


73 



where t 3 is the time at which the ball leaves the center of the catcher and 
starts rolling back off toward the edge. 

The position and the velocity of the ball as functions of time are described 
in Equations (42, 43) in the previous chapter. We also need to satisfy Equa- 
tion (50) from that chapter to insure that the ball leaves exactly backward 
(horizontally). 

Let t 4 denote the time at which the ball leaves the edge of the catcher. We 
can find the value of this time by solving the condition expressed in Equa- 
tion (50) mentioned before. This condition has been developed in Chapter 3 , 
we will then just use the result in Equation (57) (previous chapter), where we 
replace ti by £4. The value of k is adjusted so that the results are physically 
acceptable. The value of t 4 has to be greater than (t 2 + (t 4 -t 3 )). The value 
of (U - t 3 ) can be determined by solving Equation (89) at the edge where 
t = t 4 and r(t) = r x . 

Finally, ( t 3 — <,•), the transfer time between the two tracks and possibly 
some extra delay time where the ball is maintained near the rotor axis to 
insure proper orientation of the centrifuge, can now be deduced knowing 
(U - < 3 ), t 4 ) and ti. 

4.6 PRACTICAL CASE 

This is a continuation of the practical case considered previously. Let us use 
the same specifications as before for the centrifuge: 

• rx = 0.1 m, 

• ri = 595 m/s, 

• K = 5976. 
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The trajectory of the ball on the track is described by: 

r(<) = 0.1 cosh K(t - ti ) - 0.099 sinh - *,). 

Figure - 22 shows the position of the ball as a function of time. r(t) decreases 
from rj to r 0 in an interval of time (t 2 - < t ) = 0.5 ms. 

Knowing t, from the previous section, f,- = 1.3 ms, we can now deduce 
t 2 - So, the total time for the ball to reach the center of the catcher from the 
time we start spinning the centrifuges is t 2 = 1.8 ms. The position of both 
tracks is represented at this time in Figure - 23 and Figure - 24. 

We also illustrate r(t) versus t in Figure - 25. The mathematical expres- 
sion for r(t) as a function of time is: 


r(t) = 597.6 sinh K(t - *,) - 595 cosh K(t - t,). 

It starts at -r x = -595 m/s and reaches the value -4.25 m/s after 0.5 ms 
which gives a velocity value of the center of mass, according to Equation (83), 
of 104.33 m/s at r 0 near the center of the rotor ( about 6% of its initial value). 
This should help the transfer to the second track. 

Now, before we can determine how long the ball is going to be kept near 
the center, we can calculate the interval of time (t 4 - * 3 ) as explained before. 
We find: 


U — t 3 = 0.5 ms, 

still the same time as before since the ball is rolling from r 0 at rest, to the 
edge, rj. 

The next step is to calculate t 4 knowing that it has to be greater than 
2.3 ms. The result is: 


t 4 = 2.48 ms, 

with a value of k equal to 9. 
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Figure - 23. Orientation of the track on the relay at the 

instant t = t 2 = 1.8 ms. 

The ball is now at the center of the relay. 
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Figure - 24. Orientation of the 2 n< ^ track on the relay 

(acceleration part) at the instant t = t 2 = 1.8 ms. 
The ball is now being transferred to this track. 


78 




Figure - 25. Variations of r(t) from the catching 
point to the center of the relay. 
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And finally, we can deduce 1 3 : 


t 3 = 1.98 ms. 

From the instant t 2 to the instant t 3 , for 0.18 ms, the ball will be trans- 
ferred to the second track and will be kept there until the end of this time 
interval, when a new acceleration sequence begins. An illustration of the 
position of the track at this particular moment is given in Figure - 26. The 
ball is now released from the center. 

The total time for the ball to be accelerated, launched, caught, and for it 
to roll back off the catcher to its edge is 2.48 ms. It leaves with a velocity of 
1.7 km/s backward. An illustration of the rotor blade from where the ball is 
going to be launched backward at that particular time is given in Figure - 27. 

Assuming that the catcher is still at the same distance from the primary 
launcher after this first ball catching, the time of flight between the two 
centrifuges is still 0.4 ms as before. This means that the total time it takes 
for the first ball to be launched and collected back by the primary launcher 
after it accomplishes its mission is 2.88 ms. 

At this point, we can determine the momentum transferred by this first 
ball to the relay it impacted. We found here that v out = -v in , which means 
that the decrease in the ball velocity as defined in Chapter 2, becomes: 

Av = 2u ln , (90) 

and the momentum gained by the relay is: 

M{V out - V in ) = 2 mv in . (91) 

In the particular case considered, the numerical value of this momentum 
gam is 14 kg m/s, if the ball weighs 4 grams (0.5 cm radius steel ball). 
Furthermore, if the relay is 10 kg, then the velocity increment is 1.4 m/s. 
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Figure - 26. Orientation of the 2 n< * track on the relay 
at the instant t = t 3 = 1.98 ms. 

The ball is now released from the center. 
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Figure 


J 


- 27. Orientation of the 2 nd track on the relay 
at the instant t = t 4 = 2.48 ms. 

The ball now flies off the relay. 
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Now, in terms of angular momentum, we have: 


and 


Lin = — ri X rnvi, 


( 92 ) 


Lout = n X m(-fi;„). 


( 93 ) 


This means that the angular momentum of the ball before and after the 
interaction with the relay is exactly the same and A L rotor = 0. There is no 
angular momentum gained by the rotor after the launch of this first ball. This 
can be explained by the fact that the ball transfers all its angular momentum 
to the relay when it lands on the track and rolls up to the center of the rotor 
where it stops. So, in order for the ball to be accelerated on the second 
track and leave with the same velocity it impacted the relay the first time, 
the ball gets back the angular momentum it previously transferred to the 
relay. Therefore, the relay gains L in in a first step, and then loses L out in 
a second step, which is exactly the same quantity, leaving it with zero net 
angular momentum gain. This is not necessarily the case for the next balls 
launched. It depends on whether or not the ball is leaving with exactly the 
same angular momentum it had when it impacted the relay. This in turn, 
depends mainly on the position where the ball is caught on the track of the 
relay and where it exits this relay. 
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Chapter 5 


GENERALIZATION: MULTIPLE 
LAUNCH ANALYSIS 


After the detailed analysis of the launching and the catching of the first 
ball, as described in Chapter 3 and Chapter 4 respectively, we are now to 
consider multiple balls launching and catching, until our objective of acceler- 
ating the relay to the desired speed is achieved. We will start by generalizing 
the results obtained in Chapter 3 to describe the acceleration of the nth ball 
on the primary launcher. We will then describe the catching position and 
time, and the behavior of the nth ball on the relay, based on our analysis 
in Chapter 4. We will also give the velocity increments the relay is getting 
from its impact with the nth ball. All these results are to be coded in a 
computer program that would keep track of all the information conserning 
the acceleration of the spaceborne relay. 

5.1 LAUNCH OF n BALLS 

We can determine at this point the frequency of launch of the projectiles. 
The maximum frequency would be just the same as the frequency of the 
rotational movement of the centrifuge, which is in the case of the application, 
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1.59 kHz. This means that we launch one ball per revolution. However, we 
will consider launching one ball every two revolutions instead, to make sure 
that the previous ball had enough time to leave the relay, when the newly 
launched ball impacts it. When the track is in the right position to launch 
the first ball, which is t = 0.4 ms in the particular case considered before, we 
wait for exactly two revolutions (1.26 ms) and release the second ball. The 
time would then be t = 1.66 ms. The first ball would have already left the 
launcher. It will also take this second ball 0.5 ms to get to the tip of the rotor 
and it will leave at t = 2.16 ms. We can generalize this result and state that 
soon after one ball leaves the launcher, another one is released near the center 
of that launcher. The exact value of this time delay between the exit at the 
tip of the rotor of one ball and the release of the next one, is determined by 
the difference between two periods of rotation of the centrifuge and the time 
it takes the ball to reach the tip of the rotor when it is accelerated from the 
center. The n th ball will be released near the center at a time t releaae : 

t release = (0.4 -f 1.26 (u — 1)) mS, 
in the case of this application, or in general: 

^release = ^0 + 2 (n — 1) T, (94) 

and it will leave the track at 

texit ~ (0.9 + 1.26 (n — 1)) ms, 
in this particular case, which can also be generalized as: 


Uxit = <i + 2 (n - 1) T 


(95) 
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5-2 rSnsr RELAY AFTER THE LAUNCH 

Based on our kinematics calculations in Chapter 2, the velocity of the relay 
after the launch of the nth ball would be given hy Equation (8) in that 
chapter, where V, n would be the relay velocity after the (n - l)th ball leaves 
c mg back to the primary launcher. We can just put that equation in a 
loop specifying the initial velocity of the relay, which is zero (starting from 
r«t), and the program will give us the velocity of the relay after each ball 

S ' 3 smoS? N OF THE catching po- 
sition AND TIME FOR THE nTH BALL 

As w e mentioned previously, we want an initial separation of about 500 m 

t* T S V ( « Pnmary ' aUnCher “ d ‘ he relay - Theref ° rc ’ When we Equa- 
7 (7S ! (a ° P ‘ er 4) in the we keep on increasing the value of the 

:r:r: t ***** *. - such .hat the of 

a ’ t/li9ht ’ corres P° nds to a distance of at least 500 m This 
w.U set the initial distance. Then we will solve Equation (78) mentioned 

ove, to get the time at which the nth ball must be caught, and k would 
- chosen so that the result is physically acceptable. This means that the 
.me of catching needs to be at least equal to the time the nth ball exits the 
er, p us the time of flight, which can be formulated mathematically as: 

> *e xit + * flight 

For the determination of the catching position, we just solve Equation (77) 

^ exp a, ned before in Chapter 4 , and „e will find the angular position, <j, of 
s on the track at the impact point. It is important to note here that 
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one condition ban to be met when solving Equation (77). This condition is: 


1 + a tan(t), + u it/) < 0. 


(96) 


his condition was not mentioned before in Chapter 4 because it was satisfied 

‘ he “** ° f ‘ he first *»“• in “>* particular example considered. However, 
since we are generalizing the results in this chapter, we will have to check 
t is condition and make sure it is satisfied. The first ball is caught at the 
edge of the track, but, as we will see later from the results of the computer 
program, the second ball will be caught slightly away from the edge, toward 
e center. Actually, the difference in angular position, of the second ball 
catching point and the first ball is very small (one thousandth of a radian). 
This difference between two successive balls is very small, but between the 
1st and the nth, it will be quite large. The catching point will move farther 
and farther from the edge, until it becomes, eventually, near the center of 
e relay, when this one is moving at a speed comparable to the speed of 
the incoming balls. Therefore, tan te +,*) will change a great deal from 
the value calculated for the first ball, and we want to make sure that it 
a ways satisfies the condition expressed in Equation (96) above. Otherwise, 
catching cannot be achieved. The system of Equations (71,72) from the 
previous chapter, will not have a solution, and the catching condition of 
matching velocity vectors at the impact point would be impossible. It turns 
ter we run our program, that this condition (Equation (96)), limits the 
number of balls launched to 882 in the case of the specific example treated. 

is represents actually a little more than l/10th of the total number of balls 
needed to get the relay going at a velocity almost equal to the velocity of 
t e incoming balls ( v r „„ = 1.698 km/s for n = 8000, and v MI = 1.7 km/s). 
This number of balls will indeed be sufficient the accelerate the relay to more 


87 


than half the velocity of the balls. The relay velocity obtained with n = 882 
is tVe/av = 861 m/s. This is the limit of what we can achieve. A display of the 
catching position of the last ball is given in Figure - 28. We cannot possibly 
launch and catch more balls unless we change the parameters of the system, 
namely, the track shape. We tried a more circular track, with a constant 
a = 0.5, to improve the velocity of the relay, but the change would not be 
dramatic. The exit velocity of the balls would be in this case v exit = 1792 
m/s, and the maximum number of balls we can catch smoothly would be 
n = 946. The resulting relay velocity would then be u re/oy = 951 m/s. A 
display of this track is given in Figure - 29. We chose to keep the previous 
case, a = l, since the improvement is not very important compared to the 

trouble that we may have for mechanical design: machining a multiple tracks 
rotor as shown on Figure - 30. 


5.4 BEHAVIOR OF THE nTH BALL ON THE RE- 
LAY 

The nth ball will go through two different stages depending on which track 
At first, it will roll up the track from the catching position and stop 
at the center of the relay. Then, it will be transferred to the second track, 
where it will be accelerated and roll back toward the edge of the relay. It will 

fly off the track and go back to the primary launcher. We will now consider 
these two cases separately. 


5.4.1 From Catching Position to Center of Relay 

The first ball being caught at the edge of the relay, will roll for 0.5 ras to 
reach the center, which is the same as the time it takes the ball to reach the 
edge when accelerated from the center on the primary launcher. However, 
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Figure - 28. Catching position of the first and the last ball launched. 
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for the second ball, since it is caught a little farther from the edge, toward 
the center, the time needed for it to reach the center will be slightly shorter, 
and for the next ball, it will be even shorter, and so on, until the last ball 
launched. To find the exact value of this time, we need to solve Equation (88) 
in Chapter 4 , where we replace r x by r„ and -r x by r, for the general case 
of the nth ball. The equation becomes: 

r (0 = r > cosh K {t - U) + sinh K(t - *,•). (97) 

One way to solve this equation analytically would be to use the following 
identity: 

cosh 2 x — sinh 2 x = 1, 
or, 

sinh x = ±y cosh 2 x — 1. 

Now, since the argument of the hyperbolic sine in our case, K(i - *,), t being 
the time we are solving for, is positive, sinh K(t - t t ) is also positive. There- 
fore, we will chose the “+” sign over the sign in the expression above. 
Now substituting for sinh K{t - <,•) as explained above in Equation (88) men- 
tioned earlier, we get: 

r(t) = n cosh K{t ~ti) + ^r^osh 2 ^ ~t~) Z 1. ( 98 ) 

Let us denote cosh K(t - f.) by X, for simplification purposes. Then, Equa- 
tion (98) can be written as: 


r{t) = r 0 = nX + ?±y/x*-l, 
K 


or, 


ro - r t X = jVxTTT. 


(99) 


( 100 ) 
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We can square both sides of Equation (100) and rewrite it as: 


\+2r 0 r i X- + U 0. 


The solution to this equation is: 


-r 0 r t * ± y/A 

r 2 2 ’ 

K 7 r i 


where A is given by: 


r? r 2 

— — - (r 2 - r 2 ^ 
j {2 2 v t 'o; • 


We would note here that we have to chose the positive root only since X 
represents the hyperbolic cosine, which is always positive. 

The next step is to solve for the argument, K(t — tj), of the hyperbolic 
cosine. Let us call this argument x. we can write: 


e x + e~ x 


If we denote e x by Y, we can write the equation above as: 


or simply, 


* = J (y + F>’ 


K 2 - 2 XY + 1=0. 


The solutions to this equation can be written as: 


Y = X ± VX 2 - 1. 


(103) 


We have to keep in mind that Y has to be positive, since Y = e x . We 
will chose the “+” sign in the above. One could argue that the root with 
the sign is also positive. But, even when positive, it will be a relatively 


C- 1 



small number, which in turns means that the argument of the exponential, 
x = K(t - ti), would be negative. And this last statement is not physically 
acceptable. Therefore, we will chose the larger solution of the two, which is 
the one with the sign. 

We can finally write the time the needs the reach the center of the relay: 

t ~ ti = K ln [^ + VX 2 - lj , (104) 

where X is as defined in Equation (102) above. 

5.4.2 Acceleration of the nth Ball from the Center to the Edge 

Now, the ball is being accelerated from rest on the second track, where it will 
fly off the edge and travel back to the primary launcher. The orientation of 

this track is similar to the launching track on the primary launcher (a = 1). 

However, that track should be oriented such that the ball is launched back- 
ward and not forward. There is a difference of tt in the angular position 
of the ball at the exit with respect to its position on the primary launcher. 
The time it takes the ball to reach the edge is, as determined before, 0.5 ms, 
and it will be the same for all balls. The exit time can be calculated from 
Equation (57) in Chapter 3. The value of k would be chosen so that the value 
of the time obtained here is greater than the time it takes this nth ball to 
reach the center, calculated in the previous section, plus the 0.5 ms needed 
to reach the edge of the second track when accelerated from the center of the 
relay, as explained in the beginning of this paragraph. We will also note that 
we start here with a value of k = 1, instead of 0 as in the case of the primary 

launcher, to take care of the orientation of the track that has to insure a 
backward launch. 
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5.5 POSITION OF THE RELAY 


We need to determine the exact position of the moving relay when the nth 
ball ,s about to impact it. After the (n - 1 )th ball leaves the relay, this latter 
would be moving at a velocity already determined (second section of this 
chapter). Therefore, the distance it travelled from the time the (n - l)th ball 
left it, until the nth ball arrives, can be easily determined. This relay position 
should be compared to the distance that the nth ball travelled, from the time 
it exit the primary launcher, to the catching time. If the distance travelled by 
the ball is larger, as we will see from the results of the program, a mechanism 
should be included to slow down the ball for exact timing. This is one of 
the purposes of the guidance system. As we will describe it in details later 
all three directions of motion would be corrected by the aiming mechanism.’ 
Vertical adjustment will insure the catching of the ball at the desired position 
Co and a longitudinal adjustment will insure a synchronisation between the 
incoming ball and the catching time on the relay. 


5.6 ROTATIONAL ENERGY OF THE RELAY 

The last part of this kinematics and dynamics study, would be to calculate the 
change in rotational energy the relay experiences due to its impact with all 
these little projectiles. For the first ball, caught and launched from the edge, 
there is no change in rotation energy, as explained in Chapter 4 . However,' 
for the nth ball, caught at a position between the edge and the center of the 
relay, and then accelerated all the way to the edge, the relay would lose some 
rotational energy. Let us determine exactly how much energy it would lose. 
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We have: 


and 



We also have: 


and, 


' r . cos (ut { + C.) \ 
K r « + C.) / 

( T] COS (ut + Cl) | 
\ T\ sin(u>t + Ci) ) 


where t is the time the nth ball flies off the edge of the relay and all other 
quantities are as defined previously. 

Therefore, we can write the angular momentum as: 


Lin = -mr,v in sin(o)t, + C,), 

and, 

Lout = mrxVin sin(wt + Ci)- 

These two quantities can be easily calculated and the difference, L in - L outy 
which is negative, represents the decrease in rotational energy after the nth 
ball impacts and leaves the relay. We have to compile all the decreases of 
energy, from the second to the last ball, to get the total lost of rotational 
energy. We should then be able to deduce the change in angular speed of the 
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relay. We expect this change to be negligible, but we will calculate it in our 
program. 

We are now ready to write the general program that codes all the different 
steps described in this chapter. We will mention though that we added one 
more thing at the beginning of this program. It will check first the sliding 
tendency based on the roll versus skid analysis in Chapter 3. If the ball 
is expected to slide on the track, it will ask the user to change the system 
parameters. Otherwise, it will inform him that the choice of parameters 
is fine and that the ball will roll nicely on the track. It will then go on 
calculating the data needed to summarize all the derivations and visualize 
the propulsion of the spaceborne relay. 

We will include in this section a summary of the results for the first 20 
balls launched. Table - 4 will give the timing information, and Table - 5 the 
positions information. The general code written is given in the appendix. 
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Table - 4. Timing information 


^release' time at which the ball is released near the center of 
the primary launcher. 

iexit' time at which the ball flies off the launcher. 
icatch'- time at which the ball is caught on the relay. 
tbrtup- time at which the ball reaches the center of the relay. 
texrei' time at which the ball flies off the relay. 
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Table - 5. Positions information 


C.: catching position on the relay. 

Ke/oy: velocity of the relay. 

Xreiay’ position of the relay at catching time. 
■Xbaii- position of the ball at catching time. 
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appendix 


COMPUTERIZED COMMUNICATION AND CONTROL 
SYSTEM FOR SPACEBORNE CENTRIFUGAL RELAYS FOR 

SPACECRAFT PROPULSION 


One of the most important engineering problems for this project consists 
of accurately tracking the position of the projectiles and rotors and guiding each 
projectile exactly onto the appropriate guide tracks on each rotor. We have to 
control velocity vectors of the rapidly moving balls so that they encounter guide 
tracks on the spinning rotors at just the right location and time. We need to 
correct errors in lateral positions and make longitudinal adjustments to insure 
that projectiles arrive at the next relay in proper synchrony with its rotation. 

Guidance studies include the tracking mechanism (interrupt of optical 
beam) and the aiming mechanism (electromagnetic focussing including the 
electromagnetic deflection coils design and the switching circuitry). A 
schematic of the system is given in figure 1 





Figure- 1 Computerized Communication and Control System 




TRACKING MECHANISM 


Interrupt of optical beam has been identified as the tracking mechanism. 
A narrow parallel laser beam is substantially diverged by a cylindrical lens with 
extremely short focal length. The diameter of the cylindrical lens has to be 
sufficiently small, but slightly greater than the nominal width of the raw laser 
beam. The lens should be about 1 to 2cm long and 2 to 4mm in diameter. We 
will also use an array of photodiodes to detect the position of the ball. Parallel 
output arrays of the type PDA-20 will be used. The precision on the position of 
the rapidly moving ball is of the order of lO' 5 m. 

We will be using two beam-photodiode array systems to determine both 
lateral positions of the ball. Then a similar system down the track will give the 
position of the ball at a later time allowing us to determine the drift in lateral 
position and also a measure of the longitudinal speed. 

The mechanism is illustrated in figures 2 and 3. 

A second method has also been investigated although the interrupt 
beam remains the leading candidate tracking mechanism. This potential 
method consists of creating a uniform electric field that will be perturbed by the 
presence of the conducting sphere. A grid will be installed along the field and a 
measure of the potential on different points of this grid will give us the position of 
the sphere. This principle is illustrated in figure 4. 
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Figure- 2. Interrupt Beam 
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Figure- 3 . Interrupt Beam For Both Lateral Positions 



AIMING MECHANISM 


For the velocity trimming, we propose that the projectiles be directed 
through a guide barrel after release from the centrifugal accelerators to limit the 
dispersion of lateral velocities. We identified the concept of electromagnetic 
focussing as the final aiming mechanism. When the ball arrives to the 
correction circuit it will be magnetized. Its magnetic moment and the external 
magnetic field that we create will produce a potential energy. A force will then 
be acting on the sphere modifying its position. The velocity errors are expected 
to vary randomly in magnitude and in direction from one projectile to the next 
one launched. We need the switching time of the electromagnetic circuit to be 
very small in order for us to set the necessary amount of magnetic induction in 
the circuit to make the appropriate correction for each projectile. The first thing 
we need to determine is the diffusion time of the magnetic field inside the 
sphere to make sure that it will have enough time to get magnetized while 

passing through the electromagnetic deflection coil at a high velocity. For this, 
we solved Maxwell equations, 


VxH=^pJ + 

VxE+ ?f = 


1 3D 
c at 

0 


V.D = 47tp 
V.B = 0 


and got the diffusion equation, 

Jr = Vx(vx B) + ^V2 B 



The magnetic field decays away exponatially in a diffusion time: 


47tctL 2 

x = 

M-C 


For a steel ball of 1cm radius, x = 0.6ps meaning that the magnetization can be 
considered as instantaneous compared to the total time spent by the sphere on 
one set of coils (20ps) and no premagnetization is necessary. 

The magnetization of the ball can be expressed as: 


M = 


JL M. 

4n p+ 2 



where: M = magnetization 

m = magnetic moment 

V = volume of the sphere 

p = permeability of the material (Steel) 

B = magnetic induction of the external field 
and the potential energy obtained is: 


U = -m . B 


giving a force: 


where: 


F = pV y = - V.U 
p = volumic mass 
7 = accceleration 



3_ 4-1 

^~4np 4+2 


V.B 2 


This acceleration will compensate the error on the velocity: 


v = Vq - rvyt = 0 


Vo = initial error (about 0.1 m/s) 
n = number of loops 

The magnetic induction producing this acceleration is obtained when a 
current I is circulating in a loop of radius a: 


la 2 


P _ . . Icl^ 

™ 2(a 2 + z 2 ) 3/ 2 (along the z-axis) 


Using N turns instead of one in every loop will make the current needed to 

produce the same magnetic induction N times smaller. The total resistance R of 
the coils (for one loop) is given by: 


R = Ro - 

s 


where: R 0 = resistivity of the wire (copper) 

I = total length of the wire 
s = cross-section of the wire 

The voltage necessary to produce the required current is: 


V = Rl 


We have to minimize the power losses: 


P = R|2 


the inductance L of the coils, 


L 



a 2 K 


where h = height of the short solenoid 
a = radius of the winding 
K = coefficient 

and the time constant, 



We explored the tradeoff between switching speed and power losses with an 
aim to find the optimum compromise between these two parameters. An 
acceptable switching scheme has been found using the previous equations. 

For a projectile (1cm radius steel ball) launched at Ikm/s and arriving at the 
correction circuit with an error of 0.1 m/s the characteristics of the deflection coils 
are as follows. Using 10 identical sets of coils (2cm in diameter for each set), 
the magnetic induction necessary for the compensation of the velocity error is 
such that V.B 2 = 1.65 T 2 /m, the corresponding current is 10.1 A assuming 100 



turns in each set of coils. The resistance for one sef of coils is o.4n. the 
diameter of the wire is 0.5mm, Ihe inductance is O.ISmH and the time constant 
is 0.45ms. The voltage needed is 4.04V and the power losses are 40.8W 
corresponding to an increase in the temperature of 19<>K for one second. These 
studies indicate that electromagnetic deflection is feasible for the magnitude of 
velocity errors one normally obtains with a rifle barrel. 








Figure-6. Switching Circuitry 


APPENDIX 


We will present in this section the computer programs written (in c- lan- 
guage), to generate the data needed for understanding the behavior of the 
system as a function of certain parameters. We will start with RELAYVE- 
LOCITY.C, the velocity of the relay versus the number of balls launched, as 
illustrated in Figure-5 in the 2nd chapter. The second program we wrote, 
EXITVELOCITY.C, illustrates the dependency of the exit velocity of the 
balls from the primary launcher, V exit , on the track shape constant, a. An 
illustration of V ex itlVtip was given in Figure - 13. The third program written, 
ROLLSKIDW.C, represents the sliding tendency of the ball on the track as 
a function of the angular velocity uj. An illustration of the ratio A 2 /Aj, ver- 
sus u>, for four different values of the track shape constant, a, was given in 
Figure - 17. The fourth program, ROLLSKIDA.C, represents the maximized 
sliding tendency of the ball on the track as a function of the track shape 
constant a. An illustration of the ratio (A 2 /Ai) TOOa; , versus a, was given in 
Figure - 18. Finally, the general code that summarizes all the analysis done 
in Chapter 2 through Chapter 5, KIN&DYN.C, is given in this last part of 
the appendix. This code generates all the data necessary to visualize the 
concept of acceleration of a spaceborne payload. This program can be run 
for an arbitrary number of balls, but the number n was chosen equal to 20. 
The data obtained was given in Table - 4 and Table - 5. 
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/* RELAYVELOCITY . C */ 


/* THIS PROGRAM CALCULATES THE RELAY VELOCITY */ 


#include "stdio.h" 
main () 


FILE *fp; 

float V_exit=1701 ; /* this value of the exit velocity of the ball was 

calculated in the numerical application in the 
second chapter */ 

float M=10; /* mass of the relay in kg */ 

float m=4e-3; /* mass of the ball in kg */ 

float eps=m/M; 
float V[8010] ; 
int i; 


int n=8000; /* number of balls launched */ 

V[0]=0; /* relay initial velocity */ 


/* Opening of a file for the data */ 


if ( (fp=f open ("velocity . dat " , "v" ) ) —NULL) 

{ 

printf ("cannot open file/n"); 
exit(l) ; 

> 
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f* Calculation of the relay velocity */ 


for(i=0;i<n;i++) 

{ 

V [i+1] =V [i] +2*eps* (V_exit-V [i] ) ; 
fprintf (fp,'7.f\n",V[i]) ; 

> 

> 


102 



/* EXITVELOCITY . C ♦/ 


/* THIS PROGRAM ILLUSTRATES THE RELATIONSHIP BETWEEN 
THE EXIT VELOCITY OF THE BALLS FROM THE PRIMARY 
LAUNCHER, AND THE TRACK SHAPE CONSTANT. */ 


♦include "stdio.h" 
♦include <math.h> 
main() 

{ 


FILE *fp; 
float r0=le-2; 

float rl=10e-2; 
float omega=10e+3; 
float a=0.5; 
float V_tip; 
float V_exit; 

int i; 


/* distance from center of rotor to beginning 
of track (in meters) */ 

/* radius of the rotor (in meters) */ 

/* rotor angular speed (in radians per second) */ 
/* track shape constant */ 

/* rotor tip speed (in meters per second) */ 
exit velocity of the balls from the primary 
launcher */ 


float ratio [15]; 


/* Opening of a file for the data */ 

if ( (fp=f open ("exitvelocity.dat" , "v") )«NULL) 

{ 

printf ("cannot open file\n") ; 
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exit(l) ; 

} 


/* Calculation of the exit velocity of the ball */ 

V_tip=omega*rl ; 
f or(i*0; i<*10; i++) 

i 

V_exit=V_tip*sqrt ( 1+ (5 . /7) * ( 1- (rO/r 1) * (rO/rl) ) +2*sqrt (5 . /7) 
*sqrt ( 1- (rO/rl ) * (rO/rl) ) * 1/sqrt ( l+a*i*a*i) ) ; 
ratio [i] =V_exit/V_tip ; 
f printf (f p , "Xf Xf \n" , a* i , rat io [i] ) ; 

> 

> 
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/* ROLLSKIDW.C */ 


/* THIS IS AN ANALYSIS OF THE SLIDING TENDENCY 
OF THE BALL ON THE TRACK AS A FUNCTION OF 
THE ANGULAR VELOCITY OF THE ROTOR */ 


#include "stdio.h" 
♦include <math.h> 
main() 

{ 


FILE *fp; 

float rl=10e-2; /* radius of the rotor */ 

float omega; /* rotor angular speed (in radians per second) */ 

float omegaO= 1000.0; 

float a; /* track shape constant */ 

float V_tip; /* rotor tip speed (in meters per second) */ 

float KO; 

fl0at K; /+ constant appearing in the equations of motion * 

float t=0.5e-3; /* time at which ball reaches the edge 

of the track */ 

float th; 
float ratiol [25] 
float ratio2[25] 
float ratio3[25] 
float ratio4[25] 
int i; 


/* required coefficient of friction */ 
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/* Opening of a file for the data */ 


if ( (f p*f open ( "rollskidv . dat " , "w" ) ) «NULL) 

{ 

printf ("cannot open file\n") ; 
exit(l) ; 

> 


/* Calculations and writing of the data */ 


V_tip=omega*rl ; 
a*0.5; 

f or(i=l ; i<=20; i++) 

{ 

omega=omegaO*i ; 

K0*5./ (7*(l+l/(a*a))) ; 

K=omega*sqrt(KO) ; 
th=tanh(K*t ) ; 

ratiol [i] =2*a*sqrt ( l+a*a) / (5* (l+th*th)+2*sqrt (35* (l+a*a) ) *th) ; 

> 

a=l. ; 

for(i=l; i<=20; i++) 

omega=omegaO*i ; 

K0=5./ (7* (1+1/ (a*a ))) ; 

K=omega*sqrt (KO) ; 
th=tanh(K*t) ; 
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ratio2 [i] -2*a*sqrt ( l+a*a) / (5* ( l+th*th)+2*sqrt (35* (l+a*a) ) *th) ; 

> 

a=3. ; 

for(i=l;i<=20;i++) 

{ 

omega s omegaO* i ; 

K0=5 ./(7*(1+1/ (a*a) )) ; 

K=omega*sqrt (KO) ; 
th=tanh(K*t) ; 

ratio3 [i] =2*a*sqrt ( l+a*a) / (5* ( l+th*th) +2*sqrt (35* ( l+a*a) ) *th) ; 

> 

a=10 . ; 

f or (i*l ; i<=20 ; i++) 

{ 

omega=omegaO*i ; 

K0=5./(7*(l+l/(a*a))); 

K=omega*sqrt(KO) ; 
th=tanh(K*t) ; 

rat io4 [i] =2*a*sqrt ( l+a*a) / (5* ( l+th*th) +2*sqrt (35* ( l+a*a) ) *th) ; 
fprintf (fp,"Xf */.f ‘/.f */.f */.f\n", omega, ratiol [i] ,ratio2[i] , 
ratio3[i] ,ratio4[i] ) ; 

> 

> 
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/* ROLLSKIDA.C */ 


/* THIS IS AN ANALYSIS OF THE SLIDING TENDENCY 
OF THE BALL ON THE TRACK AS A FUNCTION OF 
THE TRACK SHAPE CONSTANT */ 

•include "stdio.h" 

•include <math.h> 
main() 

{ 

FILE *fp; 

float rl=10e-2; /* radius of the rotor */ 

float omega=10e+3; /* angular velocity of the rotor */ 

float a; /* track shape constant */ 

float a0=0.5; 

float V_tip=omega*rl ; 

float ratio [25]; /* required coefficient of friction */ 

int i; 

/* Opening of a file for the data */ 

if ( (f p=f open ( "rollskida . dat " , "w" ) ) *=NULL) 

{ 

printf ("cannot open file\n"); 
exit(l) ; 

> 
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/* Calculations and writing of the data */ 


for(i*0;i<*20;i++) 

{ 


a*aO*i; 

ratio [i] «2*a*sqrt ( l+a*a) /5 ; 
f printf (f p , "Xf Xf \n" , a, rat io [i] ) ; 

> 
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/* KINADYN.C */ 


/* KINEMATICS AND DYNAMICS OF THE PRIMARY LAUNCHER, 
THE RELAY, AND THE PROJECTILES. */ 

#include "stdio.h" 

# include <math.h> 
main() 

{ 

FILE *fpl; 

FILE *fp2; 


/* DEFINITION OF THE PARAMETERS OF THE SYSTEM 
(All distances are in meters and masses in kilograms) */ 


float r0=le-2; 

float rl=10e-2; 
float l=0.5e-2; 
float M=10; 
float m=4e-3; 
float omega=10e+3; 
float a»l; 
float MU_avai=0.93; 

float a_max; 


/* distance from the center of the rotor to 
the beginning of the track */ 

/* radius of the rotor */ 

/* radius of the ball */ 

/* mass of the rotor */ 

/* mass of the ball */ 

/* rotor angular speed (in radians per second) */ 
/* track shape constant */ 

/* available coefficient of friction; in this case, 
"stainless steel-stainless steel" */ 

/* maximum track shape constant allowed to 
avoid skidding */ 
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float MU_req; 

/* maximum coefficient of friction required 
to avoid skidding */ 

float epsilon=m/M; 

/* mass ratio */ 

float pi=3. 141592654; 

float T*2*pi/omega; 

/* rotor period of rotation (in seconds) */ 

float V_t ip«omega*ri ; /* rotor tip speed (in meters per second) */ 

float V_exit; 

/* exit velocity of the balls from the 
primary launcher */ 

float K0«5./(7*(l+l/(a*a))); 

float K; 

/* constant appearing in the equations 
of motion */ 

int n=20; 

/* number of balls to be launched */ 

float t_exit[25]; 

/* t_exit [n] =time at vhich the nth ball 
flies off the launcher */ 

float t_release[25] ; 

/* t_release[n]=time at which the nth ball 

is released near the center of the launcher */ 

float t_brt; 

/* time it takes the ball to roll on the 
track From rO to rl */ 

float V_relay[25] ; 

/* V_relay[n]=relay linear velocity just as 
the nth ball leaves the relay heading back 
toward the launcher */ 

float t_catch[25] ; 

/* t_catch [n] =time at which the nth ball is 
caught on the track of the relay */ 

float zeta_c[25]; 
float tg; 

/* zeta_c [n] ^angular position of the nth ball on 
the track at catching time */ 

float r_c; 

/* r_c=radial position of the ball on the track 
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float rdot.c; 
float U; 
int k; 
int i; 

float t_f light; 

float d_min; 

float d; 
float Del; 
float XI; 
float X2; 
float X; 
float t; 

float t_brtup[25] ; 


float t_exrel[25] ; 


float X_relay[25] ; 

float X_ball [25] ; 

V_relay[0]=0. ; 
K“omega*sqrt (KO) ; 


at catching time */ 

/* rdot_c“d(r_c)/dt */ 

/* rotational speed of the ball (at the tip) */ 

/* time of flight of ball between the launcher 
and the relay */ 

/* minimal initial separation between 
the launcher and the relay */ 

/* actual initial separation */ 


/* t_brtup[n]=time for the nth ball to roll up 
the track and reach the center of the relay 
when starting at the impact point */ 

/* t_exrel [n]=time at which the nth ball 
exits the relay and travels back to the 
primary launcher */ 

/* X_relay [n] “position of the relay when it 
is about to catch the nth ball */ 

/* X.ball [n] “position of the ball at 
catching time */ 

/* relay initial velocity */ 
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/* ROLL VERSUS SKID ANALYSIS */ 


a_max=sqrt ( (sqrt (l+25*MU_avai*MU_avai) -l)/2) ; 
MU_req=2*a*sqrt ( l+a*a) /5 ; 
if (MU_avai>MU_req) 

{ 


printf ("Ball will roll nicely on track. Your choice of system 
parameters is fine. You could choose 'a' as high as '/.f 
if you wanted to) . \n" ,a_max) ; 

> 

else 

{ 


printf ( Ball will skid on track. Choose another pair 
of materials to change the value of MU.avai or, 
simply choose a smaller value for 'a’ . The maximal 
acceptable value for 'a' is */.f\n",a_max) ; 
exit() ; 

} 


/* EXIT VELOCITY OF THE BALLS FROM THE PRIMARY LAUNCHER */ 
printf ("The rotor tip speed is equal to %f m/s. \n",V_tip); 

V_exit=V_tip*sqrt(l+(5./7)*(l-(rO/rl)*(rO/rl))+ 

2*sqrt(5./7)*sqrt (l-(rO/rl)*(rO/rl))*l/sqrt(l+a*a)) ; 

printf ("The speed of the ball when it flies off the track of the 
primary launcher is equal to */.f m/s. \n",V_exit) ; 
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/* FORWARD LAUNCH (PRIMARY LAUNCHER RELAY) */ 

/* Computation of the time it takes the ball to reach 
the edge of the track (distance rl) when released 
near the center at a distance rO */ 

t_brt=(l./K)*log(rl/rO+sqrt ( (rl/rO)*(rl/rO)-l)) ; 
printf ("The time it takes the ball to reach the edge of the 
track, when released near the center of the rotor 
is texit [l] -trelease [1] = y,f s . \n" ,t_brt) ; 

/* Computation of the time at which the 1st ball exits 
the launcher */ 

k*0; 

label 1 : t_exit [1] = (1 . /omega) * (atan( (-1/ a)-(V_tip/(K*sqrt ( (rl*rl) - 
(rO*rO) ))))-(l./a) *log(rl/ rO) +k*pi) ; 
while(t_exit [1] <t_brt) 

k=k+2 ; 

goto label 1 ; 

> 


/* Computation of the time at which the 1st ball 
is released near the center of the launcher */ 
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t.release [1] =t_exit [1] -t_brt ; 


/* Times at which the ith ball is released and 
exits the launcher, assuming we are launching 
one ball every 2 revolutions of the rotor */ 


for(i=l;i<=n;i++) 

{ 

t_exit[i]=t_exit[l]+2*(i-l)*T; 
t.release [i] =t .release [1] +2* (i-1 ) *T; 

> 

/* RELAY VELOCITY AFTER THE LAUNCH OF THE ith ball */ 


for(i=l;i<=n;i++) 

{ 

V.relay [ i] =V_ relay [i-1] +2* (epsilon/ ( 1+epsilon) ) * (V_exit-V relay [i-1]) • 

> 


/* DETERMINATION OF THE CATCHING POSITION AND TIME */ 
/* Catching position and time for the 1st ball */ 

U=sqrt (5 ./7)*omega*sqrt ((rl*rl)-(rO*rO)) ; 
tg* s qrt(pow( (V.exit-V.relay [0] )/U, 2 )*(l+l/pow(a, 2 ))-l) ; 

zet a.c [1] =- ( 1/a) *log ( (-U* ( 1 -a*tg) ) / ( omega*rO*sqrt ( 1+pow (a , 2) ) ) ) ; 
k-1; 
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Iabel2 : t .cat ch [1] = (1/omega) * (at an (tg) -zeta.c [1] +k*pi) ; 
while(t_catch[l] <=t_exit [1] ) 

{ 

k=k+2; 

goto label2; 

> 


/* Initial separation between the primary launcher 
and the relay */ 

t_flight=t_catch[l]-t_exit[l] ; 
d_min=V_exit*t .flight ; 

printf("The minimal initial separation is %f m. \n",d_min); 
printf ( It is better to start with a distance of at least 
500 m, which means that the ball is caught several 
revolutions later than this value of t_catch[l] 
determined above. Let us add then to t_catch[l], 
2k(pi)/omega until we reach our goal of 500 m.\n"); 

k*l; 

label3 : t .catch [1] = ( 1/omega) * (at an (tg) -zeta.c [l] +k*pi) ; 
't-fligkt=t_catch[l]-t_exit[l] ; 
d=V_exit*t_f light ; 
while(d<500) 

{ 

k=k+2 ; 

goto label3; 

> 
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printf("The initial separation is now set to %fm.\n",d); 

/* Catching position and time for the ith ball */ 

for(i=2; i<=n;i++) 

{ 

tg=sqrt (pow( (V_exit-V_relay [i-1] ) /U, 2) *(l+l/pow(a, 2) )-l) ; 
if(tg>l.) 

< 

zeta.c [i] *- ( 1/ a) *log ( (-U* ( l-a*tg) ) /(omega*rO 
*sqrt(l+pow(a,2)))); 

label4 : t.catch [i] * (1/omega) * (atan(tg) -zeta.c [i] +k*pi) ; 

while(t_catch[i]<*(t_exit[i]+t_flight)) 

{ 

k»k+2 ; 

goto label4; 

> 

> 

else 

{ 

printf ("The maximum number of balls that can be caught 
is y.d\n",i-l) ; 

exit(); 

> 

> 

/* BALL DN RELAY */ 
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/* 1st phase: from impact point to center of relay */ 


f or ( i= 1 ; i<=n ; i++ ) 

r_c=rO*exp(-a*zeta_c[i]) ; 

rdot_c= (a*r_c*omega-a*sqrt (pow (V_exit , 2) * ( l+a*a) - 
pow(a*r_c*omega,2) ))/ (l+a*a) ; 

Del=(pow((rdot_c/K) ,2))*(pow((rdot_c/K) ,2)-(r_c*r_c-r0*r0) ) ; 
Xl=((-r_c*rO)-sqrt(Del))/((pov((rdot_c/K) ,2))-(r_c*r_c)) ; 
X2=((-r_c*rO)+sqrt(Del))/((pow((rdot_c/K) ,2))-(r_c*r_c)) ; 
if (X1>0) 

if (Xl<=10) 

{ 

X*X1 ; 

> 

else 

X=X2; 

} 

> 

else 

{ 

X=X2; 

> 
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t_brtup[i] = (l . /K)*(log(X+sqrt((X*X)-l .0))) ; 

> 


/* 2nd phase: from center of relay to edge of track */ 


for(i«l;i<»n;i++) 

{ 

k=l; 

label5:t_exrel[i]=(l. /omega)* (atan((-l/a)-(V_tip/(K*sqrt((rl*rl) 

- (r0*r0) ))))-( 1 . /a) *log(rl/r0) +k*pi) ; 
t=t_catch [i] +t_brtup [i] +t_brt ; 
while (t.exrel [i] <t ) 

{ 

k*k+2 ; 

goto label5; 

> 

> 


/* POSITION OF THE RELAY */ 


X_relay[l]=d; 
f or(i=2; i<=n; i++) 

X_relay [i] =X_relay [i-1] + (V.relay [i-1] * (t_catch [i] -t.exrel [i- 1] ) ) ; 

> 


/* POSITION OF THE ITH BALL AT THE TIME OF CATCH */ 
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X_ball[l]*d; 

for(i=2;i<=n;i++) 

{ 

X.ball [i] =V_exit* (t.catch [i] -t_exit [i] ) ; 

> 

/* OPENING OF FILES FOR THE DATA */ 


if ( (fpl=f open ("time . dat" , "w") ) ==NULL) 

< 

printf ("cannot open file\n"); 
exit(l) ; 

> 

if ( (fp2=f open ("posit ion. dat" ,"w") )==NULL) 

{ 

printf ("cannot open file\n"); 
exit(l) ; 

> 


/* WRITING THE DATA */ 


f or(i=l ; i<=n ; i++) 

f printf (f pi ."/.f y.f y,f Xf %f\n",t_release[i],t_exit[i], 
t_catch[i] ,t_brtup[i] ,t_exrel[i]) ; 
fprintf (fp2,"'/,f %f y.f */,f \n" ,zeta_c[i] ,V_relay [i] , 
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X.relay [i] ,X_ball [i] ) ; 
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